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ABSTRACT 


In  this  report  we  Investigate  the  possibility  of  using  ultra-low- 
frequency  (ULF)  signals  from  airborne  loop  antennas  (I.e.,  magnetic 
dipoles)  for  air/undersea  communication.  Because  of  the  low  data  rate 
at  ULF,  communication  Is  here  understood  to  mean  the  transfer  of  short 
messages  of  high  Information  content. 

We  use  numerical  Integration  to  calculate  the  1 Hz  total  magnetic 
field  amplitudes  In  the  sea  at  depths  from  0 m to  200  m due  to  airborne 
unit  moment  vertical  and  horizontal  magnetic  dipoles  at  altitudes  from 
100  m to  10  km.  Considering  the  magnetic  moments  m attainable  with 
present  aircraft  power  and  payload  capability  (m  ■ 10^  A m*  - 10®  A m*) 
and  the  minimum  detectable  amplitude  for  a 1 Hz  signal  beneath  the  sea 
1 my),  we  conclude  that  air/undersea  communication  at  1 Hz  Is  possible 
under  the  following  Illustrative  conditions:  For  a horizontal  plane 
loop  antenna  at  3000  m altitude  and  a ULF  receiver  at  a 100  m depth, 
communication  Is  possible  for  horizontal  distances  to  10  km  for 
m • 10’  A m®  and  to  33  km  for  m * 10’  A m*.  The  corresponding  limits 
for  a vertical  plane  loop  antenna  are  13  km  to  64  km.  It  Is  also  pos- 
sible, If  desired,  to  limit  communication  to  a comparatively  small 
circular  area  directly  beneath  the  aircraft.  Sea  floor  effects  can 
alter  these  values  significantly,  particularly  If  the  receiver  Is  near 
the  floor. 
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Note:  In  this  report  we  use  ULF  (ultra-low-frequencies) 

for  frequencies  less  than  5 Hz.  Pc  1 geomagnetic 
pulsations  are  observed  In  the  upper  part  of  this 
frequency  range.  ELF  (extremely-low-frequencles) 
Is  used  to  designate  frequencies  In  the  range 
5 Hz  to  3 kHz,  VLF  (very-low-frequencles)  Is  used 
for  frequencies  In  the  range  3-30  kHz,  and  HF 
(high  frequencies)  Is  used  for  frequencies  In  the 
ranges  3-30  MHz. 
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I.  INTRODUCTION 


Submarines  have  long  played  an  important  role  In  national  defense. 

In  recent  years,  as  the  concept  of  deterrence  has  grown,  the  importance 
of  ballistic-missile  submarines  has  increased  until  they  are  now  a primary 
element  in  the  strategic-deterrent  forces  of  both  the  U.S,  and  the 
U.S.S.R.  (Scoville,  1972).  Accompanying  the  increase  in  importance  of 
submarines,  there  has  been  much  technological  innovation  to  ensure  that 
the  submarines  can  operate  for  long  intervals  at  great  depths.  Because 
sea  water  is  a conductor  and  attenuates  electromagnetic  waves,  the  in- 
creasing depth  of  operation  of  the  submarines  has  increased  the  diffi- 
culty of  communication  between  these  vessels  and  their  command  centers. 

To  help  overcome  this  difficulty,  the  frequencies  at  which  electro- 
magnetic communication  take  place,  or  are  planned  to  take  place,  have 
been  steadily  reduced. 

Unlike  ordinary  high-frequency  (HF;  frequencies  in  the  range  3 to 
30  MHz)  radio  waves,  which  do  not  penetrate  sea  water  to  any  distance, 
very-low-frequency  (VLF;  frequencies  in  the  range  3 to  30  kHz)  electro- 
magnetic waves  penetrate  a short  distance  in  sea  water  and  a number  of 
land-based  VLF  transmitting  stations  have  been  established  at  locations 
around  the  world  for  submarine  communication.  In  addition,  airborne 
VLF  systems  have  been  developed.  To  receive  the  VLF  signals,  a submarine 
operating  at  depth  can  trail  an  antenna  close  to  the  surface  and  thus 
reduce  the  attenuation  of  the  signals  produced  by  the  sea  water.  However, 
it  is  obviously  an  inconvenience,  and  it  may  also  be  difficult  or  danger- 
ous at  times,  to  deploy  these  trailing  antennas. 
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Ihere  has  been  much  research  on  the  use  of  extremely-low-frequency 
(ELF;  frequencies  In  the  range  5 Hz  to  3 kHz)  electromagnetic  waves  for 
submarine  communication.  The  Sanguine  project,  for  example,  involved 
ELF  waves  with  a carrier  frequency  of  about  45  Hz  (Wait,  1972).  Because 
of  their  lower  frequencies,  these  waves  penetrate  more  effectively  into 
sea  water  than  VLF  waves  and  they  have  other  advantages,  primarily 
propagation  related,  over  VLF  for  communication  (Bernstein  et  a1 . , 1974). 
The  data  rate  is  low  but  apparently  adequate  for  command  communication 
purposes. 

The  use  of  ultra-low-frequency  (ULF;  frequencies  less  than  5 Hz) 
electromagnetic  signals  for  undersea  communication  has  received  much 
less  attention  than  the  present  or  proposed  uses  of  VLF  and  ELF  waves. 

For  world-wide  strategic  communication,  proposed  ground-based  ULF 
antennas  suffer  from  the  same  disadvantage  as  the  ELF  antennas:  they 
have  to  be  very  large  and  require  much  power  (e.g.,  Fraser-Smith  et  a1 . , 
1972;  Grelfinger,  1972;  Davis,  Willis,  and  Althouse,  1973;  Greifinger 
and  Greifinger,  1974;  Davis  and  Willis,  1974;  Marker,  1975).  Furthermore, 
the  ULF  signals  transfer  data  at  a very  low  rate  and  they  do  not  have 
all  the  desirable  propagation  characteristics  of  ELF  waves,  Nevertheless, 
ULF  signals  penetrate  deeply  in  sea  water  with  little  loss  (the  skin  depth 
for  a 1 Hz  electromagnetic  wave  in  sea  water  is  about  250  m)  and  for  some 

purposes  they  may  provide  a useful  undersea  communication  link. 

In  this  report  we  make  a preliminary  assessment  of  the  feasibility 
of  using  ULF  signals  from  loop  antennas  (i.e.,  magnetic  dipoles)  to 
communicate  down  from  an  aircraft  to  a submerged  receiver.  This  receiver 
may  represent  a submarine.  However,  manned  submersibles,  which  are 
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being  employed  in  increasing  numbers  for  deep-sea  and  sea-bottom  re- 
search (Arnold,  1967,  Heirtzler  and  Grassle,  1976),  and  manned  instal- 
ations  on  the  sea  floor  (possibly  for  mining),  could  also  benefit  from 
a low  data  rate  communication  channel.  It  Is  also  possible  that  air- 
borne ULf^  transmitters  could  be  used  for  undersea  search  and  rescue. 

Just  as  low  frequency  electromagnetic  techniques  are  being  considered 
for  use  In  mine  rescue  (Walt,  1971). 

An  obvious  disadvantage  of  a communication  system  using  frequencies 
in  the  ULF  range  Is  the  low  rate  of  data  transfer,  particularly  In  a 
military  context  where  messages  usually  must  be  coded.  It  will  therefore 
be  Impractical  in  almost  all  cases  to  send  lengthy  textual  messages  over 
a ULF  communication  link.  However,  such  a link  should  be  adequate  for 
command  and  control  messages,  i.e.,  short  messages  of  high  Information 
content.  For  example.  It  may  be  adequate  merely  to  Instruct  a submerged 
receiver  to  rise  closer  to  the  surface  where  It  could  receive  communica- 
tions at  higher  frequencies.  It  is  for  this  specific  limited  purpose  of 
command  and  control  that  we  believe  ULF  signals  may  be  useful  for  air/ 
undersea  communication. 

If  there  Is  to  be  an  exchange  of  Information  between  the  aircraft 
and  the  submerged  receiver  It  will  probably  also  be  necessary  to  have  an 
undersea-to-air  ULF  communication  channel  In  addition  to  the  alr-to- 
undersea  channel.  However,  we  believe  the  air/undersea  channel  Is  the 
most  Important  and  also  the  most  difficult  channel  to  establish.  The 
difficulty  arises  because  of  limitations  In  the  payload  and  electrical 
power  capability  of  the  transmitting  vehicle.  A practical  ULF  trans- 
mitter Is  likely  to  be  heavy  and  require  considerable  power  and  we 
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conTidt*-  the  payload  and  power  capability  of  an  aircraft  to  be  more 
limited  than  the  same  capability  of  a submerged  ULF  transmitting  vehicle. 
Since  this  Is  basically  a feasibility  study,  we  consider  only  the  more 
limited  air/undersea  ULF  communication  channel  In  this  report. 

In  outline,  the  report  consists  of  the  following:  (a)  a chapter 

containing  a brief  discussion  of  the  theoretical  expressions  for  the  ULF 
magnetic  fields  produced  In  the  sea  by  ULF  magnetic  dipole  sources  located 
above  the  sea,  and  describing  our  numerical  Integration  technique  for 
obtaining  these  fields,  (b)  a chapter  presenting  the  results  obtained  by 
numerical  Integration  for  the  total  magnetic  field  produced  In  the  sea 
(for  depths  In  the  range  0 to  200  m)  by  unit  moment  ULF  magnetic  dipoles 
(frequency  1 Hz)  at  a range  of  elevations  above  the  sea  (elevations  In 
the  range  100  m to  10  km),  (c)  a chapter  containing  a discussion  of  the 
payload  and  power  capability  of  a P-3  Orion  aircraft,  which  we  will  take 
to  be  our  standard,  and  derivation  of  a possible  range  of  magnetic 
moments  for  an  airborne  loop  antenna,  (H)  a final  chapter  where  we  com- 
bine the  results  of  (b)  and  (c)  to  derive  the  ranges  over  which  air/ 
undersea  communication  should  be  possible  and  where  we  discuss  the 
feasibility  of  the  method  of  communication.  Because  measurements  of  the 
total  magnetic  field  can  be  omnidirectional  (unlike  measurements  of  a 
component  of  the  magnetic  field),  and  an  omnidirectional  receiving 
capability  Is  likely  to  be  a requirement  for  an  air/undersea  communication 
system,  the  calculations  and  discussion  In  this  study  are  restricted  to 
the  total  magnetic  field  produced  In  the  sea  by  elevated  ULF  magnetic 
dipole  sources. 


We  also  report  the  results  of  several  subsidiary  studies  that 
are  relevant  to  the  feasibility  study  or  to  air/undersea  cornnunl cation 
at  ULF  and  ELF  in  general.  For  example,  we  Investigate  the  relevance  of 
the  concept  of  skin  depth  at  ULF,  where  the  points  at  which  the  fields 
are  evaluated  are  often  much  less  than  a wavelength  away  from  the 
source,  and  we  show  that  the  skin  depth  still  provides  a good  Indication 
of  the  attenuation  of  ULF  fields  as  they  penetrate  sea  water.  We  also 
Investigate  the  range  of  validity  of  the  approximate  formulas  given  by 
Kralchman  (1970)  for  the  field  generated  In  the  sea  by  dipoles  above  the 
surface.  For  the  case  of  a vertical  magnetic  dipole  we  calculate  the 
effect  of  a sea  floor  on  the  fields  produced  In  the  sea  and,  finally, 
for  a very  deep  sea  and  a 1 Hz  frequency,  we  calculate  the  Impedance 
changes  produced  In  the  airborne  current  loop  by  the  presence  of  the  sea. 
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II.  METHOD  OF  CALCULATION  01  THE  UNDERSEA  FIELDS 


The  electromagnetic  field  expressions  relevant  to  the  problem  of 
communication  between  an  airborne  loop  antenna  and  a submerged  receiver 
have  been  known  for  some  time.  Their  development  can  be  traced  to  the 

pioneering  work  of  Sommerfeld  (1909).  The  expressions  involve  complex 
Integrals  (often  referred  to  as  Sommerfeld  integrals),  which  are  diffi- 
cult and  often  impossible  to  evaluate  analytically.  Thus,  until  the 
comparatively  recent  development  of  fast  and  accurate  methods  of  numer- 
ical integration  using  digital  computers,  progress  in  evaluating  the 
electromagnetic  field  expressions  for  air/undersea  communication  has 
been  severely  restricted.  An  Important  contribution  was  made  by 
Durrani  (1962,  1964),  who  derived  simplified  but  necessarily  approximate 
expressions  for  the  undersea  electromagnetic  fields  produced  by  airborne 
electric  and  magnetic  dipoles.  Additional  expressions  are  collected  in 
handbook  form  by  Kraichman  (1970).  However,  as  we  will  show,  these 
expressions  are  not  always  adequate  for  consideration  of  air/undersea 
communication  at  frequencies  in  the  ULF  range:  they  fail  at  short 

source-receiver  distances  and  are  of  uncertain  accuracy  at  intermediate 
distances.  In  this  report  we  present  the  results  of  field  calculations 
using  a numerical  integration  technique.  These  results  are  of  great 
accuracy  and  apply  at  all  frequencies  and  source-receiver  distance  at 
which  the  quasistatic  approximation  is  valid. 

In  this  chapter  we  will  first  present  the  Integral  expressions  for 
the  magnetic  fields  produced* in  the  air  by  airborne  vertical  and  hor- 
izontal magnetic  dipoles  (corresponding  to  airborne  horizontal  and 
vertical  loop  antennas,  respectively).  We  will  then  briefly  review  the 
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Figure  II-l.  Geometry  for  air/undersea  communication. 


previous  work  on  the  approximate  analytical  evaluation  of  these  expres- 
sions. Finally,  we  will  describe  our  numerical  Integration  technique 
for  evaluating  the  field  expressions. 

1 • The  Field  Expressions 

The  geometry  used  In  formulating  the  field  expressions  Is  shown  In 
Figure  II-l.  We  use  a cylindrical  coordinate  system  (p,  (p,  z),  In 
which  the  sea  (of  conductivity  <7^  ■ 4 mho/m)  occupies  the  half-space 
2^0,  and  the  complementary  region  z > o Is  free  space  (a  • o).  In 
certain  cases,  a sea  floor  (o  ■ o^)  will  be  allowed  to  occupy  the 
region  z 1 -D  otherwise  the  sea  will  be  considered  to  be  of  Infinite 
depth  (D  The  transmitter  Is  modeled  as  a harmonic  magnetic  dipole 

of  moment  me^'^'^  located  at  (h,  0,  0).  We  shall  justify  the  choice  of 
a magnetic  dipole  shortly  on  the  basis  of  the  work  of  Durrani  (1964). 

Two  dipole  orientations  will  be  treated:  a vertical  magnetic  dipole 
(VMO)  pointing  In  the  >z  direction,  and  a horizontal  magnetic  dipole 
(HMD)  oriented  along  «(i  • 0.  The  dipole  approximation  Is  valid  provided 
the  distance  R from  the  source  to  the  receiver  (located  at  p,  - d,  ti)) 
Is  much  greater  than  the  characteristic  dimension  of  the  actual  trans- 
mitter antenna.  The  expressions  for  the  VMD  were  adapted  from  the  work 
of  Ryu  et  al . (1970),  whose  fontiulatlon  permits  the  computation  of  the 
subsurface  fields  In  a layered  conducting  medium  (I.e.,  the  effects  of 
the  sea  floor  may  be  examined),  and  also  Includes  loop  antennas  of 
finite  radius  (which  permits  the  straightforward  computation  of  the 
antenna  Input  Impedance).  The  expressions  for  the  HMD  were  taken  from 
the  work  of  BaPios  (1966),  who  developed  a unified  treatment  of  the  fields 
of  ideal  dipole  sources  In  the  vicinity  of  a homogeneous  conducting  half 
space. 
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The  field  expressions  used  1n  this  work  are 


Bp(VMO)  « ^ / g^(X)  Uj  [Ug  + Uj  tanh  u^d']  (Xp)  dx  (II-1) 
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These  expressions  are  for  the  quasistatic  case,  that  is,  a11  displace- 
ment currents  are  neglected.  This  is  permissible  since,  for. all  fre- 
quencies u and  slant  ranges  R to  be  considered  here,  uR/c  « 1, 
where  c Is  the  velocity  of  light.  Equations  (1)  through  (3)  include 
a sea  bottom  of  conductivity  at  depth  D.  The  corresponding  for- 
mulas for  an  Infinitely  deep  ocean  are  obtained  by  setting  d ■ D ■ * 
or  Ug  . U3. 

2.  Review  of  Earlier  Work 

Most  previous  studies  of  air/subsurface  communication  links  have 
been  carried  out  using  asymptotic  expansions  of  the  integrals  (II-l)  - 
(II-6),  subject  to  simplifying  assumptions.  The  most  extensive  published 
analysis  Is  that  of  Durrani  (1962,  1964),  who  considers  undersea  commu- 
nication using  magnetic  dipoles,  and  derives  simplified  expressions  for 
the  electromagnetic  fields  produced  In  the  sea  (which  is  considered  to 

be  infinitely  deep)  by  airborne  vertical  and  horizontal  dipoles.  However, 
these  expressions  are  valid  only  In  the  "quasl-near"  range,  which  moves 
away  from  the  dipole  source  as  the  frequency  decreases.  Durrani  re- 
stricts his  analysis  to  frequencies  In  the  range  10  Hz  to  100  kHz,  and 
thus  does  not  consider  ULF  frequencies.  At  10  Hz,  he  calculates  the 
"quasl-near"  range  to  be  66.2  m to  4780  km,  and  estimates  his  results 
to  be  correct  to  within  0.6X  for  distances  In  the  range  562  m to  478  km, 
For  a frequency  of  1 Hz,  which  we  take  to  be  our  representative  LiLF 
frequency,  the  minimum  source-receiver  distance  R for  the  "quasi -near" 
approximation  to  be  valid  is  /T(T  times  larger  than  at  10  Hz,  i.e.,  the 
minimum  distance  for  lest  than  O.bX  error  it  about  1.8  km.  However, 
values  of  R much  lets  than  1.8  km  could  well  be  Important  for  air/ 
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undersea  communication  at  ULF  and  thus  the  simplified  expressions  used  by 
Durrani  are  not  adequate  for  this  study. 

In  the  earlier  analysis,  Durrani  (19b2)  derived  simplified  expres- 
sions for  the  electromagnetic  fields  produced  In  the  sea  by  electric 
dipoles  located  above  the  sea.  In  his  later  work  (Durrani  1964),  he 
compares  the  magnetic  dipole  results  with  those  for  the  electric  dipole 
and  concludes  that  the  magnetic  dipole  Is  more  effective  In  most  cases 
for  air/undersea  communication.  Although  his  results  are  not  directly 
applicable  to  the  ULF  range,  the  lowest  frequency  considered  by  Durrani 
Is  probably  close  enough  to  the  ULF  range  for  his  conclusion  about  the 
magnetic  dipole  to  be  valid  at  ULF.  Thus,  In  considering  the  feasibility 
of  air/undersea  communication  at  ULF,  we  have  considered  a magnetic  loop 
antenna  Instead  of  a trailing  wire  or  other  form  of  electric  dipole. 

Other  work  by  the, authors,  to  be  reported  later,  verifies  that  the 
magnetic  loop  Is  superior  to  the  electric  dipole  for  air/undersea  com- 
munication at  ULF. 

To  obtain  a complete  specification  of  the  fields  over  the  entire 
range  of  p,  h,  d,  and  w.  It  Is  necessary  to  Integrate  (II-l)  - (II-6) 
numerically.  A numerical  study  of  the  subsurface  fields  of  a horizontal 
loop  antenna  (VMD)  located  at  and  above  the  surface  of  a homogeneous 
conducting  medium  has  been  published  by  Walt  and  Spies  (1972),  who  con- 
sidered the  problem  of  communicating  with  trapped  miners.  However, 
although  they  evaluated  the  complete  field  expressions  for  a wide  range 
of  frequency  and  conductivity,  their  results  do  not  cover  the  complete 
range  of  parameters  of  Interest  In  alr/undersea  communications  at  ULF, 
and  also  their  results  are  presented  using  a form  of  normalized  variables 
which  Is  not  easily  Interp^-eted. 
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The  studies  cited  above  neglect  the  effect  of  the  ocean  floor.  As 
will  be  shown  later  In  this  report,  this  neglect  Is  justified  1f  the 
floor  Is  more  than  about  1 skin  depth  below  the  receiver,  but  the  bot- 
tom effects  become  quite  pronounced  when  the  receiver  and  the  bottom  are 
In  closer  proximity.  Although  approximate  formulas  are  given  In 
Kralchman  (1970)  for  the  fields  In  the  presence  of  the  ocean  floor,  they 
are  useful  only  at  very  great  distances  from  the  source.  To  the  authors' 
knowledge,  the  present  study  Is  the  first  to  treat  bottom  effects  In  a 
complete  fashion  and  over  the  entire  range  of  source-receiver  distances. 

In  summary,  we  have  found  that  neither  the  approximate  analytical 
approach  to  the  evaluation  of  the  field  expressions  nor  the  numerical 
Integration  approach  have  been  applied  to  the  problem  of  air/undersea 
communication  at  ULF.  Me  have  therefore  applied  our  previously-developed 
numerical  integration  technique  to  this  new  situation.  Details  of  the 
technique  are  given  In  the  following  section}  a more  complete  description 
Is  given  by  Bubenik  (1977). 

3 . Evaluation  of  the  Field  Expressions 

Two  different  techniques  were  employed  In  the  numerical  Integration 
of  the  field  expressions  (II-l)  - (II-6).  The  primary,  simpler  method 
used  Quass-Laguerre  quadrature  rules  to  reduce  the  Integrals  to  straight- 
forward summations,  as  shown  below. 

® N 

/ f(a)  e”“  da  ■ ^ w,  f(a.|) 

0 1«1 

Here,  w^  and  are  the  set  of  weights  and  abscissas  computed  with 
respect  to  the  weight  function  e"“  for  quadrature  rules  with  N ■ 2,  4, 
8,' "256.  These  quadrature  rules  were  applied  to  each  Integral  In  order 
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of  ifKredsing  N until  the  results  of  two  successive  applications 

agreed  to  within  a relative  error  limit  of  10‘“.  The  weights  and 
abscissas  were  themselves  computed  using  a fast  numerical  technique 
recently  developed  at  Stanford  University  by  Golub  and  Welsch  (1969) 
which,  unlike  previous  methods,  does  not  explicitly  require  the  roots 
of  the  Laguerre  polynomials  to  be  found. 

The  Gauss-Laguerre  technique  works  very  well  whenever  p/(h+d)  ••  2. 
Beyond  this  limit,  the  order  N required  exceeds  256  and  Increases 
rapidly  as  this  ratio  Increases,  making  this  technique  Impractical  at 
very  large  horizontal  distances.  Often,  the  approximate  formulas  given 
by  Kralchman  (1970)  are  useful  In  this  case;  but  It  sometimes  happens, 
particularly  when  the  effects  of  the  ocean  bottom  are  considered,  that 
there  exists  a regime  where  neither  the  Qauss-laguerre  technique  nor  the 
approximate  formulas  are  useful.  A second  Integration  method  was  there- 
fore developed. 

The  numerical  Integration  strategy  exploits  the  oscillatory  nature 
of  the  Integrands  In  equations  II-l  through  II-6,  which  arises  because 
of  the  presence  of  the  Bessel  functions  Jq  and  J^.  Both  the  real  and 
the  Imaginary  parts  of  the  Integrands  have  this  general  oscillatory  form. 
In  outline,  we  Integrate  In  sequence  over  the  positive  and  negative  half- 
cycles of  f(A)  until  the  contrlbiitlon  from  the  last  of  these  sub- 
integrations Is  less  than  a specified  error  bound.  In  theory,  this 
procedure  will  always  work;  however,  to  Implement  It  In  a straight- 
forward manner  Is  very  Inefficient  as  the  envelope  of  f(X)  Is  only 
weakly  damped  as  A + particularly  for  small  values  of  h and  d. 

A more  efficient  scheme  Is  to  form  a converging  sequence  of  partial  sums 
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approximating  the  integral  by  adding  the  results  of  successive  half- 
! cycle  integrations  to  each  previous  subtotal,  and  then  to  app-ly  an 

t appropriate  transformation  to  accelerate  the  convergence  of  this  sequence, 

i The  sequence  formed  by  adding  the  terms  of  alternating  signs  derived 

I 

from  the  half-cycle  sub-integrations  is  well  suited  to  the  nonlinear 

i sequence-sequence  transform  described  by  Shanks  (1955),  and  we  therefore 

i 

implemented  this  transform.  The  sub-integrals  were  computed  using 
Romberg  adaptive  quadrature.  The  Romberg  integrations  were  carried  out 
to  a relative  accuracy  of  10‘\  and  the  relative  error  bound  for  each 
complete  Integral  was  set  at  lO"**.  All  computations  were  performed  in 
double  precision  arithmetic  on  an  IBM  370/168  machine. 

The  Romberg-Shanks  procedure  was  found  to  fill  effectively  the  gap 

between  the  Gauss-Laguerre  technique  and  the  approximate  formulas.  It 
, is  useful  even  at  very  large  values  of  p/(d  + h).  However,  it  becomes 

extremely  Inefficient  at  small  values  (<1)  of  this  ratio.  The  two 
methods  are  therefore  complementary,  with  the  Romberg-Shanks  technique 
taking  over  from  Gauss-Laguerre  as  p/(d  + h)  increases  beyond  a value 
of  about  2 (Bubenik,  1977). 


II.  RESULTS  OF  CALCULATIONS 


In  this  chapter  we  present  the  magnetic  field  data  obtained  by 
numerical  integration  for  the  VMD  and  the  HMD  at  the  representative 
ULF  frequency  of  1 Hz.  The  data  were  calculated  for  unit  moment  dipoles 
and  for  a dipole  of  arbitrary  moment  the  field  values  should  be  multi- 
plied by  the  dipole  moment  to  obtain  the  corresponding  field  amplitudes. 
We  use  the  milligamma  (1  my  ■ 10"'*  Tesla)  as  our  unit  for  the  magnetic 
field.  The  purpose  of  Section  1 is  to  show  that  the  attenuation  in  sea 
water  of  the  electromagnetic  fields  from  an  elevated  harmonic  dipole  can 
be  usefully  estimated  by  using  a skin  depth  approach.  In  Section  2 and 
3 we  present  the  major  portion  of  our  calculated  magnetic  field  data 
for  the  VMD  and  HMDi  and  these  data  are  compared  in  Section  4.  Further 
results  relevant  to  the  applicability  of  Kraichman's  (1970)  approximate 
analytical  expressions,  and  the  effects  produced  by  a sea  floor,  are 
given  in  Sections  6 and  6.  Finally,  in  Section  7 we  present  data  illus- 
trating the  impedance  changes  produced  by  the  sea  in  an  airborne  current 
loop. 


1 . Verification  of  the  Skin  Depth  Approximation 

The  skin  depth  at  angular  frequency  lo  for  an  electrically  con- 
ducting but  non-magnetic  material  (permeability  Uq)  of  conductivity  o 
is  defined  by  the  equation 
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This  equation  is  typically  derived  by  considering  a plane  electromagneti 
wave  of  angular  frequency  m falling  normally  onto  the  plane  surface  of 
the  conducting  material.  It  is  found  that  the  amplitudes  of  the  fields 
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decrease  exponentially  with  depth  as  and  thus  A is  the  depth 

at  which  the  electric  and  magnetic  fields  fall  to  1/e  of  their  values 
just  inside  the  surface. 

The  skin  depth  is  a convenient  and  useful  parameter  for  character- 
izing the  attenuating  properties  of  a particular  conducting  material. 
However,  It  is  not  certain  that  the  electromagnetic  fields  from  a 
dipole  source  close  to  the  surface  of  a conducting  material  will 
exhibit  an  dependence  on  depth.  The  reason  for  the  possible 

discrepancy  is  that,  if  the  source  is  close  to  the  conducting  material 
(or  more  specifically,  less  than  a wavelength  away),  the  electromagnetic 
fields  do  not  reach  the  surface  In  the  form  of  a plane  wave.  In 
addition  to  wavefront  curvature,  the  fields  will  also  have  a signif- 
icant induction  (or  near)  field  component,  which  decays  more  rapidly 
with  distance  than  the  radiation  field  component. 

At  ULF  and  In  the  lower  part  of  the  ELF  range  the  wavelengths  of 
electromagnetic  waves  In  free  space  are  very  large.  For  example,  the 
wavelength  of  a 7,5  Hz  wave  Is  comparable  to  the  circumference  of  the 
earth,  and  at  100  Hz  the  wavelength  Is  still  large  at  30,000  km.  Thus, 
in  considering  air/undersea  communication  at  frequencies  less  than 
100  Hz,  the  source  and  receiver  are  likely  to  be  separated  by  much  less 
than  one  wavelength.  The  question  then  arises:  do  the  ULF/ELF  electro- 
magnetic fields  produced  In  sea  water  by  an  airborne  magnetic  dipole 
source  fall  off  as  below  the  sea  surface? 

To  answer  this  question,  we  calculated  the  amplitudes  of  the  total 
magnetic  field  produced  at  depths  in  the  range  0 < d < 300  m,  and  at  a 
horizontal  distance  of  p ■ 1.0  km  In  an  infinitely  deep  sea,  by  unit 
moment  horizontal  and  vertical  magnetic  dipoles  located  at  an  altituue 
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of  h = 300  m.  We  considered  the  two  frequencies  1.0  Hz  and  45.0  Hz, 
and  the  results  of  our  calculations  are  presented  in  Figures. III-l  and 
I1I-2. 

Figure  III-l  shows  the  1.0  Hz  and  45.0  Hz  total  magnetic  fields 
produced  by  the  VMD.  The  scales  In  both  this  figure  and  Figure  II 1-2 
were  chosen  to  provide  a straight  line  representation  for  an  exponen- 
tially declining  total  field.  I.e..  the  depth  scale  Is  linear,  whereas 
the  total  magnetic  field  scale  Is  logarithmic.  The  actual  calculated 
total  fields  are  given  by  the  solid  lines  which  are  straight,  Indicating 
exponential  fall  off  with  depth.  The  broken  lines  show  the  fall 

off  that  would  occur  If  the  skin  depth  concept  completely  characterized 
the  field  variation.  As  can  be  seen  In  the  Figure,  the  two  sets  of 
lines  are  close  but  do  not  coincide.  Figure  III-l  also  emphasizes  the 
great  difference  between  1 Hz  and  45  Hz  electromagnetic  signals  In 
their  ability  to  penetrate  sea  water.  The  amplitudes  of  the  1 Hz  and 
45  Hz  signals  are  nearly  the  same  Just  beneath  the  sea  surface,  but  at 
a depth  of  300  m the  amplitude  of  the  45  Hz  signal  Is  reduced  more  than 
three  orders  of  magnitude  below  the  amplitude  of  the  1 Hz  '•■Ignal. 

Figure  III-2  shows  the  1.0  Hz  and  45.0  Hz  total  magnetic  fields 
produced  by  the  HMD  for  the  two  azimuthal  angles  4)  ■ 0”  and  ■ 90®, 
The  same  conclusions  as  for  the  VMD  apply:  the  total  fields  fall  off 
exponentially  at  a rate  approximately  characterized  by  the  skin  depth 
and  there  Is  much  greater  attenuation  at  45.0  Hz  than  at  1.0  Hz. 

Although  the  results  presented  In  Figures  III-l  and  III-2  cover 
only  selected  frequencies  In  the  ULF  and  ELF  ranges,  and  a limited 
range  of  possible  receiver  depths,  source  heights,  or  source-receiver 
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TOTAL  MAGNETIC  FIELD  DUE  TO  UNIT-MOMENT  DIPOLE  (my) 

Figure  III-l.  Variation  of  the  amplitude  of  the  total  magnetic  field 
produced  at  depths  In  the  range  0 to  300  m In  the  sea  by  an  elevated 
hannonic  VMD  of  unit  moment.  The  variation  Is  shown  at.  a horizontal 
distance  of  1.0  km  from  the  axis  of  the  dipole  and  for  the  two  fre- 
quencies 1.0  Hz  and  45.0  Hz.  The  altitude  of  the  dipole  above  the 
sea  surface  Is  300  m,  and  the  sea  Is  assumed  to  be  Infinitely  deep. 

Also  shown  are  the  variations  that  would  occur  If  the  field  amplitudes 
varied  as  A-exp  (d/6),  where  A Is  the  amplitude  Just  beneath  the 
sea  surface,  d Is  the  depth,  and  6 Is  the  skin  depth. 
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TOTAL  MAGNETIC  FIELD  DUE  TO  UNIT- MOMENT  DIPOLE  (my) 

Figure  I1I-2.  Variation  of  the  amplitude  of  the  total  magnetic  field 
produced  at  depths  In  the  range  0 to  300  m in  the  sea  by  an  elevated 
harmonic  HMD  of  unit  moment.  The  variation  Is  shown  at  a horizontal 
distance  of  1.0  km  from  the  axis  of  the  dipole  and  for  the  two  fre- 
quencies 1.0  Hz  and  45.0  Hz.  The  altitude  of  the  dipole  Is  300  m 
above  the  sea  surface,  and  the  sea  Is  assumed  Infinitely  deep.  Also 
shown  are  the  variations  that  would  occur  If  the  field  amplitudes 
varied  as  A*exp  (-d/6),  where  A Is  the  amplitude  just  beneath 
the  sea  surface,  d Is  the  depth,  and  6 Is  the  skin  depth. 
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ranges,  they  nevertheless  clearly  indicate  the  validity  of  what  we  will 
call  in  this  work  the  'skin  depth  approximation':  the  ULF  and  ELF 
electromagnetic  fields  from  an  elevated  magnetic  dipole  source  are 
attenuated  In  sea  water  approximately  as  where  d Is  the 

depth  below  the  surface  and  & the  skin  depth. 

2.  Total  Magnetic  Field  Produced  In  the  Sea  by  a VHP 

Figure  I1I-3  shows  20  plots  of  the  variation  with  horizontal 
distance  of  the  amplitude  of  total  magnetic  field  produced  at  different 
depths  In  sea  water  by  an  elevated  harmonic  VMD  of  unit  moment 
(moment  »>  1 A m^)  and  a frequency  of  1 Hz,  The  sea  Is  assumed  to  be 
Infinitely  deep,  and  the  heights  h,  depths  d,  and  horizontal 
distances  p (see  Figure  II-1)  cover  ranges  which  we  believe  Include 
all  values  of  practical  Interest  at  the  present  time.  Thus,  h takes 
the  values  100,  300,  1000,  3000,  and  10,000  m;  d takes  the  values 
25,  50,  ICO,  and  200  m;  and  p varies  from  1 m to  100  km.  Note  that 
the  magnetic  field  of  the  VMD  Is  symmetrical  about  the  dipole  axis  and 
the  data  In  Figure  III-3  are  therefore  applicable  at  all  azimuthal 
angles. 

Both  the  total  magnetic  field  and  horizontal  distance  scales  are 
logarithmic  In  order  to  show  In  one  display  data  covering  many  orders 
of  magnitude.  Also,  the  logarithmic  display  converts  a variation  of 
the  magnetic  field  with  horizontal  distance  of  the  form  where 

n Is  a constant,  Into  a straight  line.  Since  we  expect  the  magnetic 
field  at  large  horizontal  distances  to  vary  as  an  Inverse  poiwer  of  the 
horizontal  distance,  the  logarithmic  display  shows  clearly  where  the 
transition  to  the  limit  of  power-law  variation  takes  place. 
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In  a later  section  we  will  compare  some  of  our  exact  calculated 
results  with  those  derived  by  using  the  approximate  expressions  of 
Kraichman  (1970).  These  expressions  are  accurate  for  large  values  of 
the  distance  R between  the  source  and  receiver.  Consequently,  in 
Figure  1II-3,  the  range  of  validity  of  Kraichman's  expressions  Is  not 
always  restricted  to  the  straight  line  sections  at  large  horizontal 
distances,  as  might  be  expected.  We  will  show  that  for  large  h (which 
automatically  implies  large  R)  the  curves  shown  in  Figure  1 1 1-3  are 
quite  closely  approximated  over  the  entire  range  of  p by  the  fields 
derived  from  the  appropriate  Kraichman  expression. 

Each  of  the  curves  in  Figure  II 1-3  can  be  considered  to  be  made  up 
of  three  sections,  depending  on  the  form  of  the  curves  and  the  horizontal 
distance:  a 'near'  section,  an  'intermediate'  section,  and  a 'distant' 
section.  In  the  near  sections  the  total  magnetic  field  does  not  change 
significantly  with  distance.  This  lack  of. variation  is  most  pronounced 
for  a dipole  altitude  of  1000  m,  where,  for  all  depths,  the  total 
magnetic  field  remains  approximately  constant  out  to  a horizontal  dis- 
tance of  approximately  1 km.  The  Intermediate  sections  represent 
regions  of  transition  in  the  distance  variation  of  the  total  magnetic 
field.  They  separate  the  near  sections  from  the  distant  sections, 
where  the  field  varies  with  distance  as  p"’’,  i.e.,  the  power  law 
variation  discussed  above.  In  Figure  I1I-3  the  curves  all  have  the 
same  slope  (n  > 4)  in  their  distant  sections,  which  agrees  with  the 
approximate  expressions  of  Kraichman  (1970).  For  small  dipole 
altitudes  (100  < h < 1000  m)  the  field  variation  in  the  intermediate 
sections  declines  monotonically  as  p increases  and  the  maximum  field 
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amplitude  is  produced  in  the  sea  directly  under  the  dipole.  However, 
for  h > 1000  m the  total  magnetic  field  has  a broad  maximum* In  the 
Intermediate  sections.  Thus,  for  a VMD  at  an  altitude  of  10,000  m,  the 
maximum  field  Is  produced  at  a horizontal  distance  of  about  b km,  and 
the  position  of  this  maximum  Is  very  nearly  Independent  of  depth  over 
the  range  shown  In  Figure  III-3, 

Comparison  of  the  curves  In  Figure  III-3  shows  that  at  large  dis- 
tances (p  > *•  10  km)  the  total  magnetic  field  produced  In  the  sea  by 
the  VMD  Increases  as  the  altitude  Increases.  Thus,  even  though  the 
field  Is  falling  off  with  horizontal  distance  as  p"**.  It  Is  possible 
for  the  field  to  Increase  as  the  range  R Increases,  For  example,  at  a 
horizontal  distance  of  20  km  and  a depth  of  200  m.  the  total  magnetic 
field  produced  by  a VMD  at  an  altitude  of  10,000  m (R  ■ 22.36  km)  Is 
greater  than  the  field  produced  by  the  same  VMD  at  an  altitude  of  100  m 
(R  ■ 20.00  km)  by  a factor  of  approximately  2b.  This  large  Increase  In 
field  occurs  even  though  the  range  Increases  by  11,8%. 

Selected  numerical  data  detailing  the  amplitudes  of  the  1 Hz  total 
magnetic  field  produced  at  a depth  of  200  m by  a unit  moment  VMD  at  two 
altitudes  (1  km  and  10  km)  are  presented  In  Table  III-l.  These  data 
will  be  used  for  comparison  with  the  results  obtained  for  the  HMD;  they 
will  also  be  used  In  the  last  chapter  to  assess  the  feasibility  of  air/ 
undersea  communication  at  frequencies  In  the  ULF  range. 

Finally,  In  Figure  111-4  we  present  data  Illustrating  the  reduction 
of  the  amplitudes  of  the  total  magnetic  field  produced  In  the  sea  as  the 
frequency  of  the  elevated  VMD  Is  Increased  above  1 Hz.  In  this  figure 
we  compare  the  variation  with  Increase  of  frequency  of  the  total  magnetic 
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TABLE  III-l 


Amplitudes  in  ml  111  gammas  of  the  1 Hz  magnetic  field  produced  at  a 
depth  of  200  m in  an  Infinitely  deep  sea  (conductivity  A mho/m)  by 
elevated  unit  moment  HMD  and  VMD  sources.  The  data  are  tabulated  for 
two  dipole  heights,  h ■ 1 km  and  10  km,  above  the  sea  surface  and  for 
selected  values  of  the  horizontal  distance  p In  the  range  0.1  to  100  km. 
Powers  of  10  are  shown  In  parentheses. 


VMO  HHD  (♦  ■ O’)  HMD  (♦«  90’) 
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TOTAL  MAGNETIC  FELD  DUE  TO  UNIT -MOMENT  DIPOLE  Imy) 


FREQUENCY  (Hz) 

Figure  III-4.  Data  Illustrating  the  effect  of  an  Increase  of  frequency 
from  1.0  to  200  Hz  on  the  total  magnetic  field  produced  on  the  sea 
surface  (d  ■ 0 m)  and  In  the  sea  (d  ■ ZOO  m)  by  an  elevated  VMD. 


field  amplitude  produced  at.  a horizontal  distance  of  1 km  both  on  the 
sea  surface  (d  = 0)  and  at  a depth  of  200  m by  a unit  moment  ,VMD  at  an 
altitude  of  300  m.  The  sea  is  once  again  assumed  to  be  Infinitely  deep 
and  the  frequency  is  varied  in  the  range  1 to  200  Hz.  For  d = 0 there 
is  only  a small  reduction  of  the  field  amplitude  as  the  frequency  is 
increased  to  200  Hz.  However,  for  d ■ 200  the  field  amplitude,  which 
is  about  half  the  field  amplitude  at  the  surface  for  a 1 Hz  frequency, 
drops  by  nearly  five  orders  of  magnitude  below  the  surface  field  ampli- 
tude as  the  frequency  is  increased  to  200  Hz. 

The  data  in  Figure  III-4  indicate  the  importance  of  the  choice  of 
frequency  for  communication  in  sea  water.  If  severe  signal  attenuation 
is  to  be  avoided,  it  appears  essential  to  choose  a frequency  such  that 
the  skin  depth  is  roughly  of  the  same  magnitude  as  the  maximum  depth  of 
the  undersea  receivers,  or  is  greater  than  this  depth.  In  Figure  III-4 
the  field  amplitude  at  200  m declines  fay  only  one  order  of  magnitude  as 
the  frequency  Increases  from  1 Hz  to  10  Hz,  i.e.,  as  the  skin  depth 
decreases  from  251  m to  80  m.  Thus,  for  receivers  operating  down  to  a 
maximum  depth  of  200  m,  a frequency  of  10  Hz  or  less  appears  most 
desirable  for  air/undersea  communication  if  signal  attenuation  in  the 
sea  is  to  be  minimized. 

3.  Total  Magnetic  Field  Produced  in  the  Sea  by  an  HMD 

Figures  III-5  and  II1-6  show  the  amplitudes  of  the  total  magnetic 
field  produced  at  different  depths  in  sea  water  for  two  azimuthal 
angles  by  an  airborne  HMD  of  unit  moment  and  a frequency  of  1 Hz.  As 
with  the  VMD,  the  sea  is  assumed  to  infinitely  deep,  and  the  heights 
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Figure  III-5.  Variation  of  the  amplitude  of  the  1 Hz  total  magnetic 
field  produced  at  depth  d 1n  an  Infinitely  deep  sea  (o  « 4 mho/m) 
by  a horizontal  magnetic  dipole  of  unit  moment  at  height  h above 
the  sea  surface.  The  fields  are  given  In  the  vertical  plane  contain 
Ing  the  axis  of  the  dipole. 
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h,  depths  d,  and  horizontal  distances  i-  also  have  the  same  values 
and  cover  the  same  ranges  as  for  the  VMD.  Unlike  the  VMD,  the  magnetic 
field  of  HMD  has  a azimuthal  angle  dependence.  The  field  amplitude  is 
the  same  for  = 0”  and  c()  « 180“,  i.e.,  directly  in  front  of  and 

directly  behind  the  dipole,  and  It  Is  also  the  same  for  <|>  = 0“  and 
(()  = 90“ . At  other  azimuthal  angles  the  total  magnetic  field  varies 
smoothly  between  these  two  extremes  (see  equations  11*4  through  II-6). 

The  format  of  Figures  III-5  and  III-6  Is  Identical  with  the  format 
for  Figure  III-3.  The  curves  can  also  be  divided  Into  'near', 
'intermediate',  and  'far'  sections,  although  their  other  properties  are 
significantly  different  from  those  of  the  VMD  curves. 

Comparing  the  HMD  curves  for  ■ 0®  and  i^i  ■ 90®  they  are  seen 
to  be  very  nearly  identical,  particularly  in  their  near  sections. 
However,  there  are  two  differences  of  note:  (1)  In  their  intermediate 
sections  the  curves  for  the  4)  - 0°  case  either  have  distinct  mlnimums 
(if  h is  large)  or  they  have  corresponding  deviations  in  slope  (if 
h is  small),  where,  In  contrast,  the  curves  for  - 90”  have  slopes 
that  vary  smoothly  with  no  mlnimums  or  deviations,  and  (2)  at  large 
distances  (p  > 20  km)  the  total  magnetic  field  amplitudes  of  the 
4)  ■ 0®  case  are  twice  as  large  as  those  for  the  4>  “ 90"  case.  For 
both  values  of  ♦,  the  curves  for  all  different  values  of  h merge 
Into  one  at  large  horizontal  distances  and  thus  the  amplitudes  of  the 
1 Hz  total  magnetic  field  become  Independent  of  the  dipole  altitude. 
Also,  the  curves  become  linear  at  these  large  distances,  with  a slope 
corresponding  to  a p~^  power  law  variation  of  the  field  amplitudes. 
This  variation  agrees  with  Kralchman's  expressions. 


qI  HMD,  <^  = 90» 

i-  h • 100  m d • 2S  m 


5 -6 
>»' 

1-8 


3 -I 

UJ 

u! 


P -4 
o 

s: 


2 3 4! 

LOG,o  HORIZONTAL  DISTANCE  p (mttars) 


Figure  1II-6.  Variation  of  the  amplitude  of  the  1 Hz  total  magnetic 
field  produced  at  depth  d In  an  Infinitely  deep  sea  (a  » 4 mho/m) 
by  a horizontal  magnetic  dipole  of  unit  moment  at  height  h above 
the  sea  surface.  The  fields  are  given  In  the  vertical  plane  con- 
taining the  dipole  but  oriented  In  a direction  perpendicular  to 
the  dipole  axis. 
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TABLE  II 1-2. 


Ratio  [B  (HMD|  <(i  ■ 0")]/[B  (VMD)]  of  the  amplitudes  of  the  1 Hz  total 
magnetic  fields  produced  at  a depth  of  200  m in  an  infinitely  deep  sea 
(conductivity  4 mho  m"‘)  by  elevated  unit  moment  dipoles.  Powers  of 
10  are  shown  in  parentheses. 
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Si'lf’ctrd  numerical  data  detailing  the  HMD  ( !■  O ')  and  HMD 

(■t  ■ 90")  magnetic  field  ampHtu  ^ for  d = ZOO  ni  and  h =.  1 km  and 
10  km  are  listed  in  Table  III-l. 

Comparison  of  the  VMD  and  HMD  Results 

The  most  Important  difference  between  the  VMD  and  the  HMD  from  the 
point  of  view  of  air/undersea  communication  at  frequencies  In  the  ULF 
range  is  the  difference  in  the  rate  of  fall -off  of  their  field  amplitudes 
at  large  distances.  For  the  VMD  the  amplitude  of  the  total  magnetic 
field  varies  as  0*“*  at  large  horizontal  distances,  whereas  for  the 
HMD  the  amplitude  varies  as  p“’.  Thus,  the  amplitude  of  the  ULF 
total  magnetic  field  produced  In  the  sea  by  an  elevated  VMD  falls  off 
more  rapidly  with  distance  than  does  the  amplitude  of  the  total  field 
produced  by  an  HMD. 

The  more  rapid  decline  of  the  VMD  fields  with  distance  Is  Illus- 
trated by  the  data  In  Table  III-Z  for  the  ratio  [B  (HMD,  <l>  = 0"]/ 

[B  (VMD)].  The  data  In  the  table  are  for  f * 1 Hz,  d -ZOO  m,  and 

h = 1 km  or  10  km.  For  p = 100  km  and  h = 1 kin  the  field  amplitude 
produced  by  the  HMD  Is  59  times  as  large  as  the  field  amplitude  pro- 
duced by  the  VMD.  The  ratio  declines  as  h Increases;  for  example, 

at  a dipole  height  of  10  km  the  ratio  drops  from  59  to  6.6. 

Although  the  VMD  field  amplitudes  fall  off  more  rapidly  with 
distance  than  the  HMD  fields,  there  are  nevertheless  several  reasons 
why  the  VMD  should  not  be  Immediately  dismissed  In  favor  of  the  HMD  for 
air/undersea  communication.  It  may  be  very  difficult,  for  example,  to 
design  airborne  loop  antennas  of  sufficiently  large  moment  to  produce 
useful  ULF  signals  at  the  large  distances  where  the  HMD  fields  are 
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',u|)erior  in  strength  (for  the  same  dipole  moment)  to  the  VMD  fields. 
Under  some  circumstances  it  may  even  be  undesirable  to  have  detectable 
dipole  fields  at  large  distances,  in  which  case  the  VMD  could  have 
advantages  over  the  HMD,  both  because  of  Its  more  rapid  fall  off  of 
field  strength  with  p and  because  of  certain  other  features  of  its 
fields  described  below.  Finally,  Installation  of  a large  current  loop 
on  an  aircraft  will  Involve  many  practical  problems  and  it  is  possible 
that  a horizontal  loop  (VMD)  will  be  easier  to  install  than  a vertical 
loop  (HMD) . 

Two  additional  features  of  the  VMD  fields  which  could  be  important 

t 

in  air/undersea  communication  are  (1)  their  lack  of  an  azimuthal  angle 
dependence  and  (2)  the  absence  of  mlnimums  such  as  occur  in  the  HMD 
(4>  = 0")  fields  (see  Figure  III-5).  The  first  of  these  features  gives 
the  VMD  some  advantage  over  the  HMD  In  that  Its  fields  would  be  more 
predictable.  However,  we  have  already  seen  that  the  HMD  fields  for 
* 0°  and  ■ 90®  are  nearly  identical  except  for  the  presence  of 
the  mlnimums  in  the  4)  = 0°  case.  These  minimums  are  an  obvious  dis- 
advantage in  the  HMD  i iilds  and  it  is  interesting  to  observe  that  they 
occur  at  points  occupied  by  maximums  in  the  equivalent  VMD  fields. 
Consequently,  the  fields  of  the  VMD  often  exceed  those  of  the  HMD  in 
the  range  of  horizontal  distances  where  the  HMD  mlnimums  occur,  and 
this  effect  produces  the  small  values  of  the  ratio  [B  (HMD,  ()i  = 0°)]/ 
[B  (VMD)]  tabulated  in  Table  III-2.  Another  situation  where  the  VMD 
fields  exceed  those  of  the  HMD  under  identical  conditions  occurs  at 
low  h (h  < 1000  m)  and  small  horizontal  distance,  where  the  VMD 
fields  may  be  nearly  twice  those  of  the  HMD. 
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Note  that  the  in1n1tnuins  occurring  in  the  HMD  (<t>  » 0°)  curves  are 
not  related  to  the  unusual  "jiggles"  observed  by  Fraser-Smith  and 
Bubenik  (1976)  In  the  ULF/ELF  magnetic  fields  produced  at  tho  sea  sur- 
face by  submerged  magnetic  dipoles.  (These  Jiggles  may  be  caused  by 
partial  cancellation  of  the  fields  propagating  in  the  "up-and-over" 
mode  by  the  fields  propagating  directly  through  the  sea  water.) 

Instead,  the  minimums  occur  close  to  the  points  where  the  fields  of  the 
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HMD  in  unbounded  free  space  are  entirely  vertical,  I.e.,  they  are  the 
field  geometry  minimums  also  observed  by  Fraser-Smith  and  Bubenik  (1976). 
Similar  minimums  would  be  observed  In  the  curves  for  the  vertical 
component  of  the  VMD  magnetic  field t they  do  not  appear  In  the  total 
magnetic  field  curves  for  the  VMD  (Figure  III-3)  because  the  horizontal 
component  is  always  dominant  oyer  the  regions  In  the  sea  where  the 
vertical  component  has  Its  minimums. 

5 . Applicability  of  Kralchman's  Approximate  Expressions 

As  we  mentioned  In  Section  III-2,  Kralchman's  (1970)  approximate 
expressions  for  the  magnetic  field  produced  In  the  sea  by  airborne 
VMD's  and  HMD's  are  accurate  provided  the  source-receiver  distance 
(I.e.,  our  range,  R)  is  large.  How  large  R should  be  is  illustrated 
by  Figure  1II-7,  which  compares  the  exact  and  approximate  1 Hz  total 
magnetic  field  amplitudes  calculated  at  a depth  of  200  m for  a VMD  at 
two  different  altitudes,  h ■ 300  m and  1500  m,  above  the  sea.  The 
exact  field  values  were  obtained  by  using  our  numerical  Integration 
technique,  and  the  approximate  values  were  obtained  by  using  the 
appropriate  Kralchman  expression. 
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It  can  be  seen  that  for  h = 1500  m there  is  very  little  differ- 


i 
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ence  between  the  curves  for  the  exact  and  approximate  field  .values, 
even  for  horizontal  distances  as  small  as  ?0  m (I.e.,  with  the  receiver 
almost  directly  beneath  the  VMD).  However,  for  h “ 300  m the  two 
curves  begin  to  differ  appreciably  once  the  horizontal  distance  becomes 
less  than  about  0.6  km.  This  suggests  that  the  fields  derived  from  the 
approximate  expressions  are  accurate  provided  the  source-receiver 
distance  Is  greater  than  about  (500*  + 600*)'/*  m,  or  781  m.  This 
result  Is  in  good  agreement  with  the  criterion  quoted  by  Fraser-Smith 
and  Bubenik  (1976)  for  the  (JLF/ELF  magnetic  fields  produced  on  the  sea 
surface  by  submerged  dipoles:  the  relative  error  In  using  Kralchman's 
approximate  expressions  Is  less  than  101»  at  a frequency  of  1 Hz  If  the 
source-recel ver  distance  Is  greater  than  750  m ('  3 6). 

6 . Effect  of  a Sea  Floor 

The  presence  of  a sea  floor  with  an  electrical  conductivity  less 
than  that  of  sea  water  will  alter  the  magnetic  fields  produced  at  points 
close  to  the  sea  floor  by  ULF  dipole  sources  above  the  sea  surface. 

The  following  simple  argument  suggests  that  the  change  will  be  an 
Increase  In  the  measured  total  magnetic  field.  It  Is  expected  that  the 
eddy  currents  Induced  in  the  sea  water  by  the  ULF  alternating  magnetic 
field  would  produce  magnetic  fields  that  oppose  the  Inducing  field  and 
thus  give  a reduced  total  field.  If  a sea  floor  Is  present  Just  below 
the  point  In  the  sea  where  the  total  field  Is  measured,  the  eddy 
currents  beneath  the  point  would  be  weaker  than  If  sea  water  alone  was 
present,  assuming  the  sea  floor  Is  less  conductive  than  sea  water. 

Thus,  the  measured  total  field  Just  above  the  sea  floor  should  be 
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TOTAL  MAGNETIC  HELD  DUE  TO  UNIT-MOMENT  DIPOLE  (m/) 


HORIZONTAL  DISTANCE  p (km) 

Figure  III-7.  Comparison  of  some  exact  magnetic  field  data  for  a VMD 
with  the  corresponding  approximate  values  derived  from  the  appropriate 
Kralchman  expression. 
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larger  than  it  would  be  at  the  same  point  if  the  sea  was  infinitely 
deep.  In  other  words,  the  ratio  B (sea  depth  = D)/B  (sea  depth  = ■) 
should  be  greater  than  one,  where  B (sea  depth  ■>  D)  is  the  total 
magnetic  field  measur  ct  the  point  above  the  sea  floor  and  B (sea 
depth  » ■«)  is  the  total  magnetic  field  measured  at  the  same  point  in 
an  infinitely  deep  sea. 

The  simple  argument  above  obviously  does  not  take  adequate  account 
of  the  phases  of  the  various  contributions  to  the  total  magnetic  field 
at  any  point  and  thus  at  large  source- receiver  distances  it  could  be 
expected  to  be  invalid.  However,  as  the  data  in  Figures  III-8  and 
II1-9  show,  the  simple  argument  gives  little  guide  to  the  actual 
variation  of  the  total  magnetic  field  in  the  presence  of  a sea  floor. 

Shown  in  the  two  figures  are  the  distance  variations  of  the  ratio 
B (sea  depth  - D)/B  (sea  depth  » »)  for  a 1 Hz  VMD  at  the  two  altitudes 
h ■ 300  m (Figure  III-8)  and  h » 3000  m (Figure  III-9),  For  each 
altitude  the  field  ratios  are  shown  for  a point  on  the  sea  surface 
(d  ■ 0 m;  bottom  panel)  and  for  a depth  of  200  m (d  » 200  m;  top 
panel).  The  sea  depth  is  varied  in  the  range  10  m to  1000  m and 
the  sea  floor  is  represented  by  a single  semi-infinite  layer  of 
conductivity  10"*  mho/m.  This  value  of  conductivity  was  chosen  for 
the  sea  floor  because  it  is  significantly  smaller  than  the  approximately 
4 mho/m  conductivity  of  sea  water  and  yet  it  is  at  least  one  order  of 
magnitude  larger  than  the  conductivi^ty  of  dry  rock,  which  is  usually 
considered  less  than  10~*  mho/m.  The  value  is  also  consistent  with 
recent  conductivity  measurements  in  the  range  10-^  to  10"*  mho/m 
obtained  by  the  Glomar  Challenger  at  depths  of  160  to  400  m in  a drill 
hole  in  the  floor  of  the  Atlantic  Ocean  (Kirkpatrick,  1977). 
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Figure  III-0.  Changes  in  the  1 Hz  total  nagnetic  field  of  a VMD  caused 
by  the  presence  of  a sea  floor  (conductivity  10'*  mho/m,  depth  D). 

The  VMD  is  taken  to  be  at  an  altitude  of  300  m and  the  fields  are 
computed  on  the  sea  surface  (d  • 0 m)  and  at  a depth  of  200  m. 
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Figure  III-9.  Changes  1n  the  1 Hz  total  magnetic  field  of  a VMD  caused 
by  the  presence  of  a sea  floor  (conductivity  10"*  mho/m,  depth  D) . 

The  VMD  Is  taken  to  be  at  an  altitude  of  3000  m and  the  fields  are 
computed  on  the  sea  surface  (d  " 0 m)  and  at  a depth  of  200  m. 
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The  data  In  the  two  figures  suggest  that  the  simple  argument  Is 
usually  valid  directly  below  the  VMD.  However,  for  horizontal  distances 
greater  than  about  400  m,  the  argument  gives  little  guidance  to  the 
actual  field  behaviour.  It  1s  Interesting  and  Important  to  observe 
that  at  large  horizontal  distances  the  VMD  field  may  be  either  reduced 
or  Increased  by  up  to  an  order  of  magnitude  by  the  presence  of  the. 
sea  floor.  The  largest  changes  appear  to  be  restricted  to  field  points 
close  to  the  sea  floor.  In  Figure  III-8,  for  example,  the  field  at  a 
horizontal  distance  of  50  km  Is  reduced  by  slightly  over  an  order  of 
magnitude  for  a point  on  the  sea  floor  at  a depth  of  200  m ( d ~ D > 

200  m) , whereas  for  the  same  point  but  for  a sea  depth  of  450  m 

(d  ■ 200  m,  0 “ 460  m)  the  field  Is  slightly  increased  by  the  presence 

of  the  sea  floor. 

These  Increases  and  decreases  of  the  total  magnetic  field  of  the 
dipole  caused  by  the  presence  of  a sea  floor  are  of  considerable  1m> 
portance  for  a1r/undersea  communication  at  frequencies  1n  the  ULF 
range  because  the  largest  changes  appear  to  occur  at  horizontal  dis- 
tances greater  than  10  km,  where  the  fields  are  becoming  small  and 
difficult  to  detect.  Thus  It  would  appear  that  the  presence  of  a sea 
floor  close  to  the  position  of  the  receiver  could  significantly  alter 
the  maximum  range  of  detection  of  the  ULF  signal  from  an  airborne 
magnetic  dipole. 

In  practice,  unless  the  undersea  ULF  receiver  1$  actually  located 
on  the  sea  floor,  the  change  In  the  maximum  range  of  detection  may  not 
be  as  marked  as  Indicated  by  the  data  In  Figures  III-Q  and  III-9.  An 
actual  sea  floor  usually  has  an  upper  layer  consisting  of  a variable 
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depth  of  unconsolidated  sediment,  which  has  a conductivity  close  to 
that  of  the  sea  water  Itself  (Kertnabon  et  a1..  1969).  For  the  receiver, 
this  sediment  layer  could  be  considered  an  inaccessible  region  with 
electrical  characteristics  similar  to  those  of  sea  water,  which  pre- 
vents the  undersea  ULF  receiver  from  coming  close  to  the  effective  sea 
floor  (which  we  suppose  has  a conductivity  that  is  much  smaller  than 
that  of  sea  water).  Thus,  the  large  changes  in  the  maximum  range  of 
detection  for  a receiver  located  on  or  close  to  the  sea  floor  are  un- 
likely to  be  observed,  or  will  not  be  as  large  as  suggested  by  the 
data  in  Figures  III-8  and  III-9. 

7 . Impedance  Changes  Produced  by  the  Sea  in  an  Airborne  Loop  Antenna 
In  addition  to  the  many  factors  already  discussed  in  this  chapter 
that  affect  the  1 Hz  total  magnetic  field  produced  in  the  sea  by  air- 
borne loop  antennas  (i.e.,  VMD's  and  HMD's),  the  mere  presence  of  the 
sea  Itself  affects  the  signal  generating  capability  of  an  airborne 
loop  antenna  by  altering  its  impedance.  If  the  loop  impedance  is  in- 
creased greatly  by  the  presence  of  the  sea,  it  may  be  difficult,  given 
the  limited  payload  capability  of  an  aircraft,  to  produce  the  1 Hz 
signal  strengths  beneath  the  sea  that  would  be  required  in  a practical 
air/undersea  communication  system.  In  this  section  vve  show  how  the 
impedance  changes  can  be  calculated  and  we  present  some  representative 
impedance  data  for  a 12  m radius  horizontal  loop  antenna  at  various 
heights  above  the  sea.  These  data  show  that  the  impedance  changes 
produced  in  the  loop  antenna  by  the  presence  of  the  sea  are  unlikely 
to  seriously  affect  the  ULF/ELF  signal  generating  capability  of  air- 
borne current  loops. 
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The  loop  input  impedance  Z (not  Including  the  ohmic  wire  resist- 
ance) is  equal  to  the  voltage  per  unit  driving  current  Induced  at  the 
terminals  due  to  the  time-varying  magnetic  flux  linking  the  loop. 

Given  a loop  of  radius  a at  a height  h carrying  a current  I of 
angular  frequency  ui,  the  loop  input  impedance  is  given  by 

Z - R + IX  - V/I  • /s  dS 

where  the  integral  of  the  vertical  magnetic  field  component  is 
evaluated  over  the  area  S of  the  loop,  R is  the  resistive  component 
of  the  impedance  and  X is  the  reactive  component.  The  integral  may 

be  rewritten  to  give  Z in  the  following  form 

Z - ioju^ira*  f (1  + J*  (Xa)  dA, 

0 

which  is  in  agreement  with  expressions  derived  by  Wait  and  Spies  (1973). 
The  first  term  in  the  above  expression  represents  the  self,  or  primary, 
impedance  of  the  loop  exclusive  of  the  ohmic  wire  resistance.  This 
self  impedance  is  Inductive.  The  second  term  represents  the  effect 
of  the  conducting  ground  (secondary  impedance)  and  contributes  both 
real  (resistive)  and  imaginary  (reactive)  parts.  As  will  be  shown 
presently,  this  secondary  Impedance  is  much  smaller  than  the  primary 
impedance,  and  it  may  thus  be  regarded  as  an  incremental  change  in 
the  total  input  impedance  of  the  loop. 

Since  there  is  no  known  general  closed-form  expression  for  the 
second  term  in  the  above  integral,  the  integration  must  be  performed 
numerically.  The  secondary  impedance  wus  evaluated  for  a 12  m loop 
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nWQUENCY  (Hl> 


Figure  IIMO.  Incremental  Increase  In  the  Input  resistance,  as  a 
function  of  frequency,  for  a single  turn  horizontal  loop  of  12  m 
radius  located  at  height  h (meters)  above  an  Infinitely  deep  sea 
of  conductivity  o - 4 mho/m. 
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at  heights  1n  the  range  0 h < 1000  m above  a conducting  "ground" 
of  conductivity  n = 4 mho/m  (I.e.,  above  the  sea)  for  frequencies 
1n  the  range  10‘*  5 f < 10^-  Hz.  The  real  and  Imaginary  parts  were 
separated  to  give  the  Incremental  loop  Input  resistance  R and  the 
reactance  X due  to  the  Induced  electric  currents  In  the  sea.  The 
reactance  was  negative;  which  we  Interpret  as  an  Incremental  decrease 
In  the  total  Input  Inductance  of  the  loop. 

In  Figures  III-IO  and  III-ll  we  present  our  numerical  results  for 
the  Incremental  loop  Input  resistance  and  the  Incremental  loop  Inductance. 
Figure  III-IO  shows  the  Incremental  Increase  In  the  Input  resistance, 
as  a function  of  frequency,  for  a single  turn  horizontal  loop  of  12  m 
radius  located  at  height  h above  an  Infinitely  deep  sea  of  conductivity 
a ■ 4 mho/m.  Figure  III-ll  shows  the  absolute  value  of  the  Incremental 
decrease  In  the  Input  Inductance,  as  a function  of  frequency,  for  the 
same  12  m radius  loop.  As  would  he  expected,  for  a given  h,  both 
the  Incremental  resistance  and  the  absolute  value  of  the  Inductance 
Increase  rapidly  with  frequency,  and  for  a given  frequency,  these 
quantities  decrease  as  the  height  of  the  loop  Is  Increased.  Even  for 
a loop  on  the  sea  surface  ( h - 0 m),  the  values  of  Incremental 
resistance  and  Inductance  are  small  and,  as  we  will  now  show,  they 
are  also  likely  to  be  very  small  compared  to  the  actual  ohmic  resistance 
and  self  Inductance  of  any  practical  loop  that  Is  used  for  air/undersea 
communication. 

Suppose  the  12  m radius  loop  consists  of  a single  turn  of  AWG  10 
copper  wire.  The  ohmic  resistance  of  this  loop  Is  0.25  n and  Its 
self  Inductance  L,  calculated  from  the  formula  L ■ u^a  [Hn  (^)  - J-], 
where  b Is  the  wire  radius.  Is  0.143  mH.  Although  these  values  are 
small,  the  Incremental  Input  resistances  and  Inductances  caused  by  the 
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Figure  III«11.  The  absolute  value  of  the  incremental  decrease  In  the 
Input  Inductance,  as  a function  of  frequency,  for  a single-turn 
horizontal  loop  of  12  m radius  located  at  height  h above  an 
Infinitely  deep  sea  of  conductlvUy  a ■ 4 mho/m. 


sea  water  are  many  orders  of  magnitude  smaller.  For  example,  if  the 
12  m loop  is  being  driven  at  a 1 Hz  frequency  and  is  at  an  altitude 
of  1000  m above  the  sea  surface,  the  Incremental  loop  input  resistance 
is  close  to  6 X 10‘^’  n,  which  is  clearly  negligible  compared  to  the 
0.25  ohmic  resistance  of  the  loop.  For  the  same  loop,  the  incremental 
loop  inductance  caused  by  the  sea  is  about  -3  x 10"’  mH,  which  also  is 
negligible  compared  to  the  0.143  mH  self-inductance  of  the  loop. 

The  numerical  data  plotted  in  Figures  III-IO  and  IIl-ll  are  for  a 
single-turn  loop  and  should  be  scaled  in  proportion  to  the  square  of 
the  number  of  turns  for  a multi  turn  loop.  The  self-inductance  of  an 
isolated  multiturn  loop  also  scales  as  the  square  of  the  number  of 
turns  and  thus,  tur  a multi  turn  loop  above  the  sea,  the  incremental 
loop  inductance  will  always  remain  small  compared  to  the  self-inductance 
even  if  the  number  of  turns  is  very  great.  The  ohmic  resistance  of  a 
multi-turn  loop,  on  the  other  hand,  scales  as  the  number  of  turns,  and 
it  would  be  possible  for  the  incremehtal  input  resistance  to  approach 
and  even  exceed  the  ohmic  resistance  for  a multiturn  loop  with  a great 
many  turns  above  the  sea.  In  the  present  case,  however,  *the  number  of 
turns  required  for  the  resistances  to  become  comparable  (-10*’)  is 
so  great  that  the  loop  would  impractical.  Thus,  we  conclude  that 
impedance  changes  produced  by  the  sea  in  an  airborne  loop  antenna  are 
unlikely  to  be  significant  for  air/undersea  communication. 
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IV.  POSSIBLF  MAGNETIC  MOMENT^  OF  AIRBORNE  ULF  CURRENT  LOOPS 


1 . Introduction 

The  amplitude  of  the  undersea  electromagnetic  field  generated  by  an 
airborne  dipole  antenna  Is  proportional  to  the  moment  of  the  dipole. 

For  a magnetic  dipole  formed  by  a conducting  loop  enclosing  an  area  A, 
the  moment  Is  given  by  m * lA,  where  I Is  the  total  current  through 
a cross  section  of  the  loop.  Thus  the  moment,  and  therefore  the  maximum 
range  of  communication,  of  a practical  'oop  antenna  is  limited  Ly  the 
physical  size  of  the  antenna  and  by  the  total  ’Available  oower.  These 
constraints  are  particularly  rigorous  on  an  aircraft,  whore  the  size, 
weight,  and  available  power  are  subject  to  definite  restrictions.  It 
Is  the  purpose  of  this  chapter  to  estimate  the  maximum  practical  dipole 
moment  attainable  using,  as  a representative  standard  aircraft,  the 
LocKheed  P-3  Orion.  We  shall  consider  both  conventional  and  superconduct- 
ing loops. 

2 . Constraints  Imposed  by  the  Aircraft 

According  to  Jane's  All  the  World's  Aircraft,  the  Lockheed  P-3  Orion 
is  the  current  American  standard  aircraft  for  locating  and  tracking  sub- 
marines. Essentially  a modified  Electra  airliner,  it  can  carry  a 9100  kg 
expendable  load,  Is  very  maneuverable,  and  can  operate  at  altitudes  rang- 
ing from  100  m to  8600  m.  Its  wingspan  is  30.4  m and  Its  overall  length 
Is  35.6  m.  The  electrical  power  supply  Includes  thr'^.*  60  kVA  400  Hz  AC 
generators  and  an  additional  60  kVA  auxiliary  power  unit. 
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3 . Loop  j\ntennd  Design  Consi d orations 

We  begin  with  the  expression  for  the  loop  magnetic  moment, 

m =•  AI 
* AN1 

where  i is  the  current  in  each  of  N conductors  passing  through  a 
cross  section  of  the  loop.  The  loop  area  A is  limited  by  the  size  of 
'i  e aircraft  and  may  be  regarded  as  cor':t,ant.  We  therefore  design  around 
the  Ni  product,  which  should  be  as  large  as  possible  within  the  con- 
straiii!*!  of  allowable  conductor  weight  and  generator  capacity. 

The  power  fed  to  the  loop  tern.inals  is  a complex  quantity,  repre- 
senting the  real  power  dissipated  in  the  loop  conductors  and  the  reactive 
power  due  to  the  loop  inductance.  Since  the  loop  inductance  is  very 
difficult  to  determine  (its  computation  is  complicated  by  the  presence 
of  the  metal  aircraft,  and  it  seems  best  determined  experimentally  from 
a model),  we  estimate  a reasonable  power  factor  to  be  0.5.  Thus,  given 
a conductor  of  DC  resistance  p per  unit  length,  the  total  (real  plus 

reactive)  power  drawn  from  the  generator  is  2i*  paN/ff  volt-amperes, 
where  u = 2/f~  for  a circular  loop  and  u = 4 for  a square  loop. 

The  Ni  product  can  be  apportioned  in  several  ways:  (1)  A single 

conductor  wound  N times  around  the  loop.  This  configuration  gives  a 
relatively  high  inductance,  requires  a high  driver  voltage  but  low 
terminal  current,  and  is  subject  to  complete  failure  if  the  conductor 
breaks.  (2)  N single  conductors  wound  about  the  loop  and  connected 
in  parallel  at  the  terminals.  This  results  in  a lower  inductance,  low 
driving  voltage  but  high  terminal  current,  and  is  relatively  unaffected 
by  breaks  in  a limited  number  of  conductors.  (3)  N^  single-conductor 
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loops  of  (N/N.j)  turns  connected  in  parallel.  This  compromise  gives  an 
intermediate  inductance  and  allows  the  terminal  voltage  and  current  to 
be  adjusted  to  the  generator's  characteristics  by  varying  . All  of 
the  above  combinations  dissipate  the  same  real  power,  and  have  the  same 
weight  for  a given  loop  moment  and  conductor  size.  For  reference,  the 
weight  of  the  conductors  in  a loop  antenna  is  W oNv^,  where  W is 
the  weight  per  unit  length  of  the  conductor. 

The  conductor  current  i is  constrained  by  power  supply  considera- 
tions, and  also  by  the  conductor  size.  Electrical  codes  and  practice 
fix  an  upper  limit  on  the  order  of  0.1  I^,  where  is  the  wire 
fusing  current  given  by  • Kd’''*,  where  K is  a constant  corres- 
ponding to  the  wire  material,  d is  the  wire  diameter  given  by 
d ■ 0,0825  exp  (-  0.116  AWG)  m,  and  AWQ  is  the  American  Wire  Gauge  size 
of  the  conductor  (ITT,  1968,  p.  4-55). 

Considering  the  large  number  of  variables  involved,  there  appears 
to  be  no  single  design  formula  leading  to  an  optimum  loop  of  greatest 
moment  given  the  area,  weight,  and  power  limitations.  A computer  optimiz- 
ation approach  would  probably  be  profitable  in  the  design  of  an  actual 
antenna,  but  this  was  not  attempted  in  this  study.  Instead,  a digital 
computer  was  used  to  determine  the  characteristics  of  many  possible 
designs,  subject  to  the  appropriate  constraints. 

For  the  P-3  Orion  we  take:  (1)  Total  conductor  weight  ■ expendable 

load  ■ 9100  kgi  (2)  Power  available  • 120  kVA  (2  generators  in  service, 

1 in  reserve) i real  power  available  ■ 120/2  ■ 60  kWj  (3a)  Loop  area 
■ 452  m*,  corresponding  to  a loop  strung  from  nose  towingtip  to  tail 
to  wingtip  to  nose  (VMO);  and  (3b)  Loop  area  ■ 9.3  m*,  corresponding 
to  a loop  wound  around  the  fuselage  (HMD). 
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The  estimated  attainable  loop  moments  under  these  conditions  are 
listed  in  Table  IV-1.  It  1s  apparent  that  the  maxiiruiri  attainable  moment 
of  a non-suoerconductlng  loop  probably  lies  In  the  range  10'  - 10®  A m'. 


Table  IV-1.  The  estimated  maximum  feasible  moments  of  ULF  trans- 
mitting loops  on  a P-3  Orion  aircraft,  using  non-superconductIng 
loops  and  a continuous  square  wave  drive  current.  Conductor 
current  1 ^ 0.03  I#;  smaller  currents  give  lesser  moments, 
while  larger  currents  require  excessive  power. 


COPPER  WIRE 


M (A  m') 

N 

1(A) 

AWG 

Weight  (kg) 

Power  (kW) 

VMD 

1.0  X 10' 

2262 

10.0 

10 

8.997 

62 

HMD 

1.7  X 10“ 

17752 

10.0 

10 

9,000 

63 

ALUMINUM  WIRE 

M (A  m') 

N 

1(A) 

AWG 

Weight  (kg) 

Power  (kW) 

VMD, 

1 .4  X 10' 

738 

42.1 

0 

8,996 

58 

HMD 

2.3  X 10® 

5792 

42.1 

0 

9,000 

58 

Using  superconducting  technology,  the  real  power  required  to  drive 
the  loop  is  almost  negllglblej  and  the  limitation  Is  the  reactive  power 
which  the  generator  can  supply  to  the  essentially  Inductive  load. 
However,  In  this  case  the  horizontal  loop  size  probably  must  be  reduced 
to  allow  for  the  necessary  cryogenic  cooling  apparatus,  thus  reducing 
the  loop  area  A and  the  vertical  dipole  moment  for  a given  N1  product. 
Nevertheless,  the  attainable  moment  may  be  quite  large;  an  existing 
superconducting  loop  of  1.5  ni  diameter  (Morrison  and  Dolan,  1973)  gives 
a moment  of  10®  A m'  at  45  Hz  using  a conservative  current  drive,  and 
weighs  less  than  1 ton.  Extrapolating  this  data  to  higher  currents  and 
larger  sizes.  It  seems  possible  that  a superconducting  transmitter  loop 
could  attain  moments  of  10®  - 10®  A m*. 


4 • Conclusions 

Vertical  magnetic  dipole  moments  on  the  order  of  10’'  - 10“  A 
appear  to  be  attainable  using  non-superconductIng  loops  on  a P-3  air- 
craft with  a continuous  square  wave  drive  current.  The  practical 
horizontal  dipole  moment  Is  nearly  an  order  of  magnitude  less.  These 
moments  could  probably  be  more  than  doubled  If  the  antenna  were  driven 
by  a pulsed  current  waveform  at  a duty  cycle  of  50*  or  less.  A super- 
conducting loop  Is  expected  to  attain  a moment  of  10‘  - 10“  A m’',  owing 
to  Its  reduced  power  requirements.  The  gain  offered  by  superconducting 
technology  would  be  especially  pronounced  In  the  case  of  the  horizontal 
dipole  (loop  In  vertical  plane),  whose  area  would  not  need  to  be  sub- 
stantially reduced  In  order  to  accommodate  the  cryogenic  cooling 
apparatus. 
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4.  Conclusions 


Vertical  magnetic  dipole  moments  on  the  order  of  10^  - 10®  A m^ 
appear  to  be  attainable  using  non-superconductIng  loops  on  a P-3  air- 
craft with  a continuous  square  wave  drive  current.  The  practical 
horizontal  dipole  moment  Is  nearly  an  order  of  magnitude  less.  These 
moments  could  probably  be  more  than  doubled  If  the  antenna  were  driven 
by  a pulsed  current  waveform  at  a duty  cycle  of  50%  or  less.  A super- 
conducting loop  Is  expected  to  attain  a moment  of  10"  - 10®  A m^,  owing 
to  Its  reduced  power  requirements.  The  gain  offered  by  superconducting 
technology  would  be  especially  pronounced  In  the  case  of  the  horizontal 
dipole  (loop  In  vertical  plane),  whose  area  would  not  need  to  be  sub- 
stantially reduced  In  order  to  accommodate  the  cryogenic  cooling 
apparatus. 
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V.  FEASIBILITY  OF  AIR/UNDERSEA  COMMUNICATION 
WITH  AIRBORNE  LOOP  ANTENNAS 


By  coinb1n-‘ng  the  results  obtained  In  Chapter  III  (results  of  cal- 
culations of  the  undersea  fields)  and  Chapter  IV  (magnetic  moments  of 
possible  airborne  current  loops)  we  are  now  able  to  estimate  the  horizon- 
tal ranges  over  which  air/undersea  communication  at  1 Hz  with  airborne 
loop  antennas  may  be  feasible. 


1 . Minimum  Detectable  1 Hz  Magnetic  Field  Amplitude 

Conventional  ULF  recording  systems  typically  use  Induction  loop 
antennas  as  sensors.  These  Induction  loops  may  consist  of  long  solenoids 
with  a core  of  ferromagnetic  material  (e.g.  mumetal)  or  large-diameter 
air-cored  colls.  An  Induction  loop  antenna  measures  a single  component 
of  the  magnetic  field  and  measurement  of  the  total  magnetic  field 
requires  three  of  these  antennas  In  a mutually  orthogonal  configuration. 

The  Internal  noise  level  in  the  Induction  loop  systems  varies  from 
system  to  system.  A representative  noise  level  Is  0.3  mr//fTz  in  the 
0.1  to  10  Hz  band  (Buxton  and  Fraser-Smith,  1974).  This:  noise  can  be 
reduced  substantially  by  using  superconducting  magnetometers  as  ULF 
receivers  (Buxton  and  Fraser-Smith,  1974;  Fraser-Sm1th  and  Buxton, 

1975).  However,  while  there  may  be  some  advantage  to  reducing  the  In- 
ternal noise  below  0.3  mY/vffif,  substantial  further  reduction  Is  probably 
not  Justified  In  the  presence  of  the  natural  ULF  noise  background.  The 
level  of  the  natural  ULF  noise  background  Is  greatly  variable:  It  varies 

In  particular  with  time  of  day,  frequency,  location  on  the  earth's 
surface,  and  with  the  state  of  geomagnetic  disturbance.  The  middle 
latitude  study  by  Fraser-Smith  and  Buxton  (1975)  shows  the  nighttime 


background  level  in  the  0.1  to  10  Hz  band  varying  from  about  20  inf/ZHz' 

at  0.1  Hz  to  about  0.4  at  10  Hz  during  geomagnetical  1y  quiet 

periods  (the  background  varies  approximately  as  f‘^  throughout  the 
frequency  range  10"'*  to  100  Hz).  During  the  day,  and  during  intervals 
of  geomagnetic  disturbance,  the  background  noise  level  increases.  At 
most  times  there  is  a minimum  of  activity  in  the  3 to  7 Hz  band,  and  the 
noise  level  in  this  band  Is  typically  less  than  1.0  n)Y//Ffz.  At  1 Hz 
the  noise  level  is  usually  less  than  2.0  my/vlTz  (the  noise  level  can 
be  up  to  100  times  larger  during  the  course  of  relatively  infrequent 
Pc  1 geomagnetic  pulsation  events  (Jacobs,  1970)).  Because  the  noiso 
level  increases  rapidly  as  the  frequency  drops  below  0.1  Hz,  it  appears 
that  the  best  frequency  band  for  ULF  air/undersea  communication  is  1 to 
6 Hz. 

For  ULF  measurements  in  the  sea,  and  particularly  near  the  sea 
surface,  the  motion  of  sea  waves  in  the  geomagnetic  field  provides  an 
additional  source  of  ULF  noise  (e.g.,  Warburton  and  Caminiti,  1964, 
Rodney,  1975).  The  amplitude  of  both  this  noise  and  of  the  natural  ULF 
background  declines  with  depth  in  the  sea,  although  the  rate  of  decline 
of  the  wave-induced  noise  may  not  be  the  same  as  the  attenuation 

that  is  to  be  expected  for  the  natural  background  (c.f.  Section  III-1). 
Measurements  of  the  wave-induced  ULF  noise  do  not  appear  to  have  been 
made  in  the  1 to  6 Hz  frequency  band.  However,  measurements  at  lower 
frequencies  suggest  that  the  noise  occurs  predominantly  at  frequencies 
below  0.1  Hz,  and  that  at  depths  of  100  m or  greater  the  wave- induced 
noise  in  the  1 to  6 Hz  bend  has  an  amplitude  that  is  much  less  than  1 my 
(Maclure  et  al . , 1964) . 
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Thus,  for  submerged  ULF  receivers  at  rest,  we  conclude  that  a 1 in^ 
total  magnetic  field  amplitude  from  an  airborne  loop  antenna  .operating 
at  a 1 Hz  frequency  can  probably  be  detected  with  little  difficulty  In 
the  presence  of  both  Internal  ULF  receiver  noise  and  natural  ULF  back- 
ground noise,  provided  the  receivers  are  sufficiently  deep  In  the  sea 
to  avoid  wave-induced  noise  at  the  operating  frequency.  Large  amplitude 
Pc  1 geomagnetic  pulsation  events  may,  on  comparatively  rare  occasions, 
produce  a total  magnetic  field  amplitude  much  greater  than  1 my  at  1 Hz. 
However,  the  characteristics  of  the  signal  from  the  airborne  loop  antenna 
win  be  so  different  from  those  of  the  Pc  1 pulsation  event  (e.g.  Fraser- 
Cmlth,  1977)  that  there  should  be  little  difficulty  distinguishing  the 
two  ULF  signals. 

If  the  submerged  receiver  Is  In  motion,  the  receiver  sensitivity 
may  be  decreased  substantially  by  a variety  of  motion-induced  noises. 

ULF  receivers  using  solenoids  or  colls  as  sensors  (Including  superconduc- 
ting magnetometers)  ore  particularly  susceptible  to  motion  Induced  noise. 
Because  they  respond  with  high  sensitivity  to  components  of  the  magnetic 
field,  small  angular  displacements  of  the  sensors  caused  by  vibration  or 
other  motion  can  Introduce  large  noise  components.  One  possible  solution 
to  this  problem  Is  to  use  total  field  magnetometers  as  ULF  receivers. 
However,  present  total  field  magnetometers  do  not  achieve  the  high 
sensitivity  of  the  component  magnetometers:  they  are,  In  particular, 
at  least  several  orders  of  magnitude  less  sensitive  than  present  super- 
conducting magnetometers. 

There  have  been  two  Important  recent  developments  In  response  to 
the  motion  noise  problem:  (1)  the  Naval  Research  Laboratory  (NRL)  Is 
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investigating  the  possibility  of  combining  operation  of  a three-axis 
superconducting  magnetometer  on  a towed  stabilized  platform  with  a 
signal  processing  method  that  sums  the  squares  of  the  three  component 
measurements  to  obtain  a total  field  output  (Wolfe  et  al.,  1974;  Davis 
et  al.,  1977),  and  (2)  magnetic  field  loop  antennas  In  the  form  of  long 
towed  cables,  with  superior  noise  characteristics  for  this  type  of 
antenna,  have  been  designed  and  tested  (Burrows,  1975).  These  develop- 
ments were  stimulated  by  the  proposal  to  use  ELF  signals  for  submarine 
communications,  I.e.,  by  the  proposed  Sanguine  and  Seafarer  ELF  communi- 
cation systems,  and  the  frequencies  of  Interest  lie  In  the  approximate 
range  30  to  130  Hz.  These  frequencies  are  substantially  higher  than  the 
ULF  frequencies  (~  1 Hz)  being  considered  for  air/undersea  communication 
In  this  report,  and  thus  1t  Is  not  always  possible  to  relate  the  progress 
being  made  In  motion-noise  reduction  to  the  ULF  range.  Nevertheless,  NRL 
has  demonstrated  considerable  progress  toward  achieving  their  required 
receiver  sensitivity  of  10'^  mY//Hz  In  the  30  to  130  band.  This  sen- 
sitivity Is  over  an  order  of  magnitude  greater  than  the  sensitivity  of 
the  best  conventional  ULr  receivers  at  rest.  Considering  the 
flat  frequency  response  of  superconducting  magnetometers,  It  Is  reason- 
able to  assume  that  the  same  high  sensitivity  could  be  achieved  In  the 
0.1  Hz  to  10  Hz  band. 

We  therefore  conclude  that  a conventional  ULF  receiving  system  at 
rest  In  the  sea  can  probably  easily  detect  a 1 my  total  magnetic  field 
amplitude  at  a frequency  of  1 Hz  provided  wave-induced  noise  can  be 
Ignored.  For  a ULF  receiving  system  In  motion,  motion-induced  noise  can 
become  an  Important  problem  If  the  ULF  receiving  system  uses  Induction 
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loop  antennas  as  sensors.  However » recent  progress  at  NRL  suggests 
that  a towed  three-ax1s  superconducting  magnetometer  can  be  designed  to 
measure  a 1 my  total  magnetic  field  amplitude  at  1 Hz.  Thus,  1n  tne  re- 
mainder of  this  chapter  we  will  adopt  a 1 my  total  magnetic  field  ampli- 
tude that  can  be  detected  by  a submerged  ULF  receiver. 

2 . Maximum  Distances  for  Alr/Undarsea  Communication  at  1 Hz 

Based  on  the  discussion  In  Chapter  4,  we  adopt  10^  to  10^  A m^ 
as  a possible  range  of  maximum  moments  for  airborne  VMD's  and  HMD's. 

Given  a minimum  detectable  total  magnetic  field  amplitude  on  the  order 
of  1 my  (at  a 1 Hz  frequency),  as  discussed  In  the  previous  section,  we 
can  now  calculate  the  ranges  of  horizontal  distances  over  which  air/ 
undersea  communication  at  1 Hz  using  airborne  current  loops  will  be 
possible. 

For  most  airborne  loop/undersea  receiver  configurations,  the  range 
of  horizontal  distances  can  be  specified  completely  by  a single  value  of 
the  distance,  p max,  which  Is  the  maximum  horizontal  distance  at  which 
the  airborne  loop  produces  a 1 my  total  magnetic  field  amplitude  (for  a 
1 Hz  frequency)  at  the  submerged  ULF  receiver.  Mith  the  exception  of  a 
few  high-altitude  VMD  configurations,  the  total  magnetic  field  amplitude 
at  any  horizontal  distance  less  than  p max  will  be  greater  than  1 my. 

In  the  exceptional  VMD  cases,  the  total  magnetic  field  amplitude 
falls  below  1 my  for  some  p min  < p max.  As  shown  by  the  data  In  Figure 
III-3,  a high-altitude  VMD  (3  km  ^ h <:  10  km)  may  produce  a detectable 
total  magnetic  field  amplitude  over  the  more  limited  range  p min  < p £ p 
max,  where  p min  > 0.  For  example,  suppose  a O.b  x 10*  A m^  VMD  Is 
located  at  an  altitude  of  10,000  m and  Its  1 Hz  total  magnetic  field 
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Is  being  measured  by  a ULF  receiver  at  a depth  of  200  la.  Scaling 
distances  approximately  from  the  appropriate  curve  In  Figure 'I II-3, 
total  magnetic  field  amplitudes  greater  than  1 my  are  produced  only  In 
the  annular  area  beneath  the  VMD  defined  by  800  m < p < 16  km.  In  this 
case  p min  - 800  m and  there  Is  a "dead"  area  directly  beneath  the  VMD 
where  the  total  magnetic  field  amplitude  Is  slightly  less  than  the 
detectable  level  of  1 my. 

It  Is  clear  from  Figure  III-3  that  such  "dead"  areas  beneath  the 
transmitter  occur  only  under  limited  conditions,  Therefore,  unless 
specifically  stated  otherwise.  In  the  follow  discussion  It  Is  to  be 
assumed  that  a receivable  field  Intensity  exists  at  all  values  of 
p < p max. 

Referring  again  to  Figure  1 1 1-3  for  the  VMD,  we  assume  our  ULF 
receiver  Is  located  at  a depth  of  200  m (maximum  attenuation),  and  we 
find  the  maximum  horizontal  distances  for  total  magnetic  fields  of 
10"’  my  and  10"’  my,  which  correspond  to  a 1 my  field  for  VMD  moments 
of  10’  A m^  and  10’  A m’ , respectively.  The  range  of  p max  thus  obtained 
for  the  VMD  at  an  altitude  of  3000  m Is  9.0  km  to  30.2  km.  (These 
values,  and  the  following  values  for  p max,  are  derived  from  the  numer- 
ical data  used  to  prepare  Figures  III-3  and  III-5;  they  are  not  scaled 
directly  from  the  figures).  For  a VMD  at  an  altitude  of  10,000  m the 
total  magnetic  field  produced  by  a 10’  A m’  moment  Is  everywhere  less 
than  1 my  and  Is  therefore  judged  to  be  unreceivable.  However,  If  the 
moment  of  the  dipole  Is  10’  A m’,  the  total  magnetic  field  exceeds  1 my 
out  to  a horizontal  distance  of  39.0  km. 
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These  results  show  the  Importance  of  having  as  large  a magnetic 
moment  as  possible  for  the  VMD  if  large  values  of  p max  are  desired. 

Also,  the  results  show  that,  for  a given  VMD  moment,  p max  can  be  In- 
creased (with  a limit  set  by  the  dipole  moment)  by  Increasing  the  alti- 
tude of  the  VMD. 

If  we  now  assume  our  ULF  receiver  Is  located  at  a depth  of  100  m, 
the  horizontal  distances  over  which  air/undersea  communication  at  1 Hz 
Is  possible  are  Increased  over  the  corresponding  distances  for  a ULF 
receiver  at  a depth  of  200  m.  However,  because  of  the  comparatively 
low  attenuation  of  1 Hz  electromagnetic  fields  In  sea  water,  the  In- 
crease In  distance  achieved  by  halving  the  receiver  depth  Is  not  as  sub- 
stantial as  may  be  expected.  For  VMD's  moments  1n  the  range 
10'  A m*  to  10*  A m^  located  at  an  altitude  of  1000  m,  the  range  of 
maximum  distances  for  air/undersea  communication  Is  8.0  km  to  26.0  km, 

The  corresponding  range 'of,  p max  for  the  VMD  at  an  altitude  of  3000  m 
Is  10.0  km  to  33.1  km.  When  the  VMD  Is  at  an  altitude  of  10,000  m,  a 
10^'  A m'  moment  Is  Insufficient  to  produce  detectable  total  magnetic 
field  amplitudes  (p  max  ■ 0).  However,  when  the  VMD  moment  Is  10*  A m^, 
detectable  total  magnetic  field  amplitudes  are  produced  at  all  horizontal 
distances  out  to  p max  ■ 43.2  km.  Comparing  the  maximum  distances  for 
communication  using  a 10’  A m*  VMD,  we  see  that  the  p max  of  39.0  km  for 
a ULF  receiver  at  a depth  of  200  m Is  Increased  to  43.2  km  If  the  ULF 
receiver  Is  at  100  m:  an  Increase  of  only  11  percent. 

As  noted  In  Section  III-4,  the  total  magnetic  field  produced  In  the 
sea  by  an  elevated  VMD  falls  off  more  rapidly  at  large  distances  than 
does  the  field  produced  by  an  HMD.  Me  would  therefore  expect  the 
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maxiiiiuiii  horizontal  distances  obtained  for  air/undcrsea  communication 

with  an  HMD  to  be  larger  than  those  obtained  for  a VMD  of  the  same 
moment.  In  fact,  under  the  same  conditions  as  for  the  VMD  above,  the 
values  of  (>  max  are  approximately  doubled  for  an  HMD  source  if  the  total 
magnetic  fields  are  measured  In  the  t|)  » 0°  plane.  For  example,  suppose 
the  ULF  receiver  Is  at  a depth  of  100  m and  Is  located  directly  ahead 
of  a 10^  A m^  HMD  (l.e.,  4)  ■ O'").  The  maximum  horizontal  distance  for 
an  air/undersea  communication  at  1 Hz  In  this  example  Is  64.0  km,  and 
this  value  of  f>  max  applies  for  the  HMD  at  any  altitude  In  the  approx- 
mate  Interval  1 km  to  10  km.  For  comparison,  we  obtained  p max  - 26.0 
km  for  a 10®  A m^  VMD  at  an  altitude  of  1 km  and  p max  * 43.2  km  fo]"  the 
same  VMD  at  an  altitude  of  10  km  (In  both  cases  the  ULF  receiver  was  at 
a depth  of  100  m).  For  the  same  HMD  (<li  ■ 0“)  source,  but  with  the  ULF 
receiver  at  a 200  m depth,  p max  ■ 56.2  km  for  all  HMD  altitudes  In  the 
approximate  Interval  1 km  to  10  km. 

These  data  emphasize  a further  distinction  between  the  VMD  and  HMD 
fields  that  1s  not  discussed  In  Section  III-4,  l.e.,  at  large  distances 
the  HMD  total  magnetic  field  does  not  Increase  with  Increase  of  the 
dipole  altitude,  whereas,  as  we  have  Just  observed.  Increasing  the 
altitude  of  the  VMD  can  Increase  the  total  magnetic  field  It  produces  at 
large  distances.  In  Figures  I1I-5  and  III-6  the  lack  of  an  altitude 
dependence  of  the  HMD  fields  at  large  distances  appears  as  a merging  of 
the  curves  at  large  p. 
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3 . Limited  Ranges  for  A1 r/Undersea  Communi cation  at  1 Hz 

In  most  communication  systems  the  priina^'y  object  is  to  obtain  the 
largest  possible  range  of  communication  for  a given  Input  power.  However, 
under  some  circumstances  it  could  be  desirable  to  limit  the  range  of 
communication.  We  have  already  noted,  in  the  previous  section,  how  air/ 
undersea  communication  at  1 Hz  with  a VMD  could  be  restricted  to  an 
annular  area  beneath  the  VMD.  However,  the  conditions  for  producing  this 
limited  range  of  communication  Involved  very  specific  altitudes  and 
moments  for  the  dipole  source. 

It  Is  also  possible,  as  the  data  In  Figure  1II-3,  III-5,  and  III-6 
show,  to  limit  communication  to  a small  circular  a'c.*  -'.Irectly  beneath 
the  VMD  or  HMD  source.  Because  the  conditions  for  producing  this 
limited  circular  area  of  communication  are  much  less  strict  than  those 
for  producing  the  limited  annular  area  of  communication,  It  Is  conceiv- 
able that  the  ability  to  restrict  air/undersea  communication  to  a small 
area  could  have  practical  use. 

The  first  requirement  for  obtaining  the  small  area  of  communication 
is  for  the  aircraft  to  fly  at  a low  altitude.  For  example,  suppose  the 
ULF  source  is  a VMD  and  that  its  altitude  is  300  m.  At  this  low  altitude 
the  flat  portions  of  the  curves  In  Figure  III-3  extend  out  only  to 
horizontal  distances  on  the  order  of  200  to  400  m.  Now,  if  the  dipole 
moment  is  reduced  (this  is  the  second  requirement)  from  its  assumed 
value  in  the  range  10’’  A m^  - 10’  Am'  to  a value  near  10’  A m',  the 
data  In  Figure  III-3  show  that  air/undersea  communication  can  be  limited 
to  the  flat  portions  of  the  h 300  m curves,  i.e.,  communication  is  re- 
stricted to  a circular  area  of  radius  approximately  400  m directly 


63 


beneath  the  VMD.  This  result  is  essentially  independent  of  the  depth 
of  the  ULF  receiver  in  the  range  0 m to  200  m,  although  the  1 Hz  signal 
will  be  close  to  1 my  and  only  just  observable  at  a depth  of  200  m for 
a 10®  A 111^  VMD.  Examination  of  Figures  III-5  and  III-6  shows  that  the 
same  limited  area  of  communication  can  be  obtained  with  an  HMD  source. 

Limitation  of  air/undersea  communication  to  a circular  area  of 
only  400  m radius  beneath  an  airborne  dipole  source  may  not  be  desirable 
in  practice  In  view  of  the  additional  requirements  Imposed  on  the  air/ 
undersea  communication  link  by  the  speed  of  the  aircraft  and  the  low 
data  rate  for  communication  at  1 Hz.  The  above  example  Is  only  illus- 
trative, The  Important  conclusion  remains,  however,  that  air/undersea 
communication  at  1 Hz  using  airborne  loop  antennas  can  be  restricted  to 
submerged  ULF  receivers  In  a circular  area  of  known  size  directly  be- 
neath the  aircraft,  and  that  the  size  of  this  circular  area  can  be  con- 
trolled within  wide  limits  by  variation  of  (a)  the  dipole  moment  and 
(b)  the  altitude  of  the  aircraft. 

4.  Conclusions 

In  this  report  we  have  presented  arguments  and  data  to  justify 
two  basic  conclusions: 

(1)  Horizontal  or  vertical  airborne  loop  antennas 
with  magnetic  moments  In  the  range  10^  A m^ 
to  10’  A m^  appear  to  be  attainable  with 
present  aircraft  power  and  payload  capability. 

It  Is  assumed  that  the  current  in  these  air- 
borne loop  airborne  loop  antennas  can  be 
switched  or  alternated  at  a 1 Hz  frequency. 
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(i1)  An  alternating  total  magnetic  field  with  a 
frequency  of  1 Hz  and  an  amplitude  of  1 my 
can  be  detected  beneath  the  sea  by  using 
present  induction  loop  or  superconducting  ULF 
receivers  at  rest  or.  with  some  further  develop 
ment.  by  using  a towed  stabilized  three-axis 
superconducting  magnetometer  system. 

Following  these  two  basic  conclusions,  we  further  conclude: 

(iii)  Air/undersea  communication  at  1 Hz  is  possible 
using  horizontal  or  vertical  airborne  loop 
antennas  under  the  following  Illustrative 
conditions.  First,  for  a ULF  receiver  at  a 
100  m depth  in  the  sea  and  a VMO  (horizontal 
loop  antenna)  at  an  altitude  of  3000  m above 
the  sea,  air/undersea  communication  at  1 Hz 
is  possible  for  horizontal  distances  in  the 
range  0 to  10  km  for  a VMD  moment  of  10’  A m* 
or  in  the  range  0 to  33  km  tor  a VMD  moment 
of  10’  A in^ 

Second,  for  an  HMD  (vertical  loop  antenna) 
at  an  altitude  of  3000  m and  a ULF  receiver 
located  at  a depth  of  100  m directly  ahead 
or  behind  of  the  HMD,  a1r/undersea  communi- 
cation at  1 Hz  is  possible  for  horizontal 
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distances  in  the  range  0 to  13  km  for  an 
HMD  moment  of  lO’'  A m“  or  in  the  range 
0 to  64  km  for  an  HMD  moment  of  10^  A m'. 

Note  that  the  maximum  distance  of  communi- 
cation for  a VMD  can  he  Increased  (a)  by 
Increasing  the  height  of  the  dipole  (within 
a limit  set  by  the  size  of  the  dipole 
moment) and  (b)  by  reducing  the  depth  of  the 
ULF  receiver,  whereas  at  a given  azimuthal 
angle  the  maximum  distance  of  communication 
for  an  HMD  can  only  generally  by  increased 
by  reducing  the  depth  of  the  ULF  receiver. 

(1v)  For  both  the  VMD  and  the  HMD  it  Is  possible,, 
to  limit  air/undersea  communication  at  1 Hz 
to  a comparatively  small  circular  area 
directly  beneath  the  dipole  source.  The 
radius  of  this  circular  area  can  be  adjusted 
within  wide  limits  by  varying  (a)  the  d1po.le 
moment  and  (b)  the  altitude  of  the  dipole, 
and  It  can  be  as  small  as  400  m for  a 10^  A m^ 
VMD  at  an  altitude  of  300  m (for  a ULF 
receiver  at  any  depth  In  the  range  0 m to 
200  m) . 

For  the  VMD,  It  Is  also  possible  to  limit 
alr/undersea  communication  at  1 Hz  to  an 
annular  area  centered  on  the  point  directly 
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below  the  dipole  source.  In  this  case  a 
circular  "dead"  area  Is  produced  directly 
below  the  dipole  where  the  magnetic  field 
is  too  weak  to  be  measured.  The  conditions 
for  producing  the  limited  annular  area  re- 
quire dipole  altitudes  in  the  approximate 
range  1 km  to  10  km  and.  for  a given  dipole 
altitude,  involve  a highly  restricted  range 
of  dipole  moments.  Because  of  these  re- 
strictions, this  mode  of  air/undersea  com- 
munication is  considered  less  likely  to  have 
a practical  use  than  the  mode  involving  com- 
munication to  a limited  circular  area  directly 
beneath  the  dipole. 

(v)  The  maximum  range  for  air/undersea  communi- 
cation at  1 Hz  is  smaller  for  a VMD  than  it 
is  for  an  HMD  of  the  same  moment.  However, 
the  VMD  total  magnetic  field  does  not  have 
either  (a)  an  azimuthal  angle  dependence  or 
(b)  field  geometry  miniriiums.  In  a situation 
where  range  is  not  a primary  consideration, 
the  more  regular  behaviour  of  the  VMD  field 
may  give  the  VMD  an  advantage  over  the  HMD 
for  air/undersea  communication. 
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(vi)  rior  ULF  receivers  on  or  close  to  a sea  floor 

low  electrical  conductivity  (o  £ 10“®  mho/ni)i 
our  calculations  indicate  that  the  total 
mal^gnetic  field  produced  by  an  airborne  loop 
anjbenna  at  horizontal  distances  greater  than 
abbut  10  km  can  be  increased  or  decreased  by 
up  I to  an  order  of  magnitude  by  the  presence 
ofjthe  sea  floor  (compared  with  the  total 
madnetic  field  that  would  be  produced  at  the 
sane  points  in  an  infinitely  deep  sea).  We 
thiirefore  believe  that  sea  floor  effects 
reiuire  further  study. 


(vii)  Chji 

an 

ant 


nges  in  the  impedance  of  airborne  loop 
ennas  caused  by  the  proximity  of  the 
ennas  to  the  sea  water  are  likely  to  be 
pletely  negligible. 


com 

Throughout  this  report  we  have  used  1 Hz  as  a representative  ULF 
frequency.  Consideration  of  the  natural  ULF  noise  background,  the  low 
data  rate  at  iIhz;  and  the  low  attenuation  of  1 to  6 Hz  electromagnetic 
fields  in  a 200  m maximum  depth  of  sea  water  lead  us  to  the  following 
conclusion; 


(viii)  The  more  favorable  ULF  frequencies  for  air/ 
undersea  communication  lie  In  the  range  1 to 
6 Hz. 
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Finally,  in  sunmary,  we  conclude: 


(1x)  Airborne  loop  antennas  operating  at  frequencies 
In  the  ULF  range  and  with  maximum  moments  In  the 
range  10’  A m’  to  10®  A m*  offer  a feasible  means 
of  transmitting  command  and  control  messages 
(I.e.,  short  standardized  messages  of  high  In- 
formation content)  to  submerged  ULF  receivers  at 
depths  In  the  range  0 to  200  m and  at  horizontal 
distances  In  the  approximate  range  0 to  60  km. 
Possible  Improvement  In  the  sensitivity  of  the 
submerged  ULF  receivers,  particularly  submerged 
ULF  receivers  In  motion,  and  the  provision  of 
greater  power  for  the  loop  antennas  could  con- 
ceivably Increase  the  maximum  horizontal  range 
to  about  100  km. 

In  a situation  where  secure  communication  Is 
a primary  requirement  (and  not  maximum  hor- 
izontal range),  ULF  communications  from  an 
airborne  loop  antenna  can  be  limited  to  sub- 
merged ULF  receivers  In  a restricted  area 
directly  beneath  the  loop  antenna.  The  size 
this  area  can  be  controlled  within  wide 
limits  by  variation  of  the  altitude  and 
moment  of  the  loop  antenna. 
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b.  Recommendations 


The  preceding  conclusions  are  based  largely  on  the  results  of  our 
calculations  of  the  1 Hz  total  magnetic  field  produced  In  the  sea  by 
airborne  loop  antennas  and  on  Information  that  1s  sometimes  incomplete 
or  not  In  an  entirely  suitable  form  for  use  In  the  present  study. 
Further  progress  In  air/undersea  communication  at  ULF  will  require 
additional  field  calculations  and  new  Information  as  suggested  In  the 
following  recommendations. 

(I)  Design  of  a compact,  lightweight,  superconducting 
current  loop  with  a high  magnetic  moment 

(m  i 10’  A m’)  that  can  be  switched  or  alter- 
nated at  frequencies  In  the  ULF  range  would 
represent  a major  advance  for  air/undersea 
communication  at  ULF.  Further  studies  of  high- 
moment  superconducting  current  loops  Is 
strongly  recommended. 

(II)  Continued  development  of  highly-sensitive  towed 
total  field  magnetometers,  particularly  for  use 
at  ULF,  would  complement  the  effort  to  develop 

an  airborne  loop  antenna  of  high  magnetic  moment. 
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(111)  Measurements  of  the  ULF  background  magnetic 
noise  In  the  sea  at  different  depths  (In  the 
range  0 to  200  m)  for  various  sea  states  and 
for  different  levels  of  natural  geomagnetic 
activity  are  desirable-  Measurements  both 
on  and  above  the  sea  floor  are  needed,  as  well 
as  measurements  In  a very  deep  sea  (no  sea 
floor  effect). 

(1v)  The  research  recommended  In  (111)  would  be 

strengthened  by  a parallel  theoretical  effort 
to  derive  spectrums  of  the  ULF/ELF  electro- 
magnetic fields  generated  at  different  depths 
In  the  sea  for  representative  sea  states. 

(v)  Additional  calculations  of  the  effects  pro- 
duced by  a realistic  sea  floor  on  the  ULF 
fields  generated  In  the  sea  by  elevated 
dipole  sources  would  help  establish  the 
Importance  of  a sea  floor  In  air/undersea 
communication. 
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